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ABSTRACT 

We use a sample of rest-frame UV selected and spectroscopically observed galaxies at redshifts 
1.9 < z < 3.4, combined with ground-based spectroscopic Ha and Spitzer MIPS 24 /jm data, to derive 
the most robust measurements of the rest-frame UV, Ha, and infrared (IR) luminosity functions (LFs) 
at these redshifts. Our sample is by far the largest of its kind, with over 2000 spectroscopic redshifts 
in the range 1.9 < z < 3.4 and ~ 15000 photometric candidates in 29 independent fields covering a 
total area of almost a square degree. Our method for computing the LFs takes into account a number 
of systematic effects, including photometric scatter, Lya line perturbations to the observed optical 
colors of galaxies, and contaminants. Taking into account the latter, we find no evidence for an excess 
of UV-bright galaxies over what was inferred in early z ~ 3 LBG studies. The UV LF appears to 
undergo little evolution between z ~ 4 and z ~ 2. Corrected for extinction, the UV luminosity density 
(LD) at z ~ 2 is at least as large as the value at z ~ 3 and a factor of ~ 9 larger than the value at 
z ~ 6, primarily reflecting an increase in the number density of bright galaxies between z ~ 6 and 
z ~ 2. Our analysis yields the first constraints anchored by extensive spectroscopy on the infrared 
and bolometric LFs for faint and moderately luminous (Lboi ^$ 10 12 L Q ) galaxies. Adding the IR 
to the emergent UV luminosity, incorporating independent measurements of the LD from ULIRGs, 
and assuming realistic dust attenuation values for UV-faint galaxies, indicates that galaxies with 
Lbol < 10 12 Lq account for « 80% of the bolometric LD and SFRD at z ~ 2 — 3. This suggests 
that previous estimates of the faint-end of the Lboi LF may have underestimated the steepness of the 
faint-end slope at Lboi < 10 12 L . Our multi- wavelength constraints on the global SFRD indicate 
that approximately one-third of the present-day stellar mass density was formed in sub-ultraluminous 
galaxies between redshifts z = 1.9 — 3.4. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high redshift — galaxies: 
luminosity function — galaxies: starburst — infrared: galaxies 
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galaxies is a central component 
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Shaplev et all 120011 : iGiavalisco et al.lll996l) . While this 
period in the Universe was perhaps the most active in 
terms of galaxy evolution and accretion activity, it was 
not until recently that advances in detector sensitivity 
and efficiency, the increased resolution and light- 
gathering capability afforded by larger 8-10 meter class 
telescopes, and a number of new powerful imagers and 
spectrographs on space-based missions such as HST, 
Spitzer, and Chandra, allowed for the study of large 
numbers of galaxies at z ~ 2. These developments have 
prompted a spate of multi-wavelength surveys of high 
redshift galaxies from the far-IR/submm to IR, near-IR, 
optical, and UV, enabling us to examine the SEDs of 
star forming galaxies over much of the 7 decades of 
frequency over which stars emit their light e ither directly 
or in di rectly through dust process ing (e.g., ISteidel et"al] 



2003 1 [2001 IDaddi et all l2004blfal: iFranx et al 
van Dokkum et alJ 120031 120041: kbraham e~ 
Chapman et al.ll2005l: lSmai1l2003h . 

The first surveys that efficiently amassed large samples 
of high redshift galaxies used the observed U n GlZ colors 
of galaxies to identify those with a deficit of Lyman con- 
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tinuum flux (e.g.. lSteidel et al.ll995h in the U a band (i.e., 
U "drop-outs") for galaxies at z ~ 3. Those initial results 
have been a dapted to select galaxies at higher redshifts 
(z > 4; e.g.. IBouwens et alj|2005l 12004 iDickinson et all 
12004 iBunker et all 12004 lYan et all l2003h _ang^oder- 
ate redshifts (1.4 < z < 3; lAdelberger et all 12004 
iSteidel et alj [2004) . Combining th ese high redshift re - 
sults with those from GALEX (e.g. JWvder et al.ll2005l ). 
we now have an unprecedented view of the rest-frame 
UV properties of galaxies from the epoch of reionization 
to the present, perhaps the only wavelength for which 
star forming galaxies have been studied across more than 
~ 93% of the age of the Universe. The accessibility of 
rest-frame UV wavelengths over almost the entire age of 
the Universe makes rest-frame UV luminosity functions 
(LFs) useful tools in assessing the cosmic star formation 
history in a consistent manner. 

The foray of observations into the epoch around z ~ 2 
has occurred relatively recently, and with it have come 
various determinations of the UV LFs at t hese epochs 
(jGabasch et all 12004 iLe Fevre etldl |2005[ ). Unfortu- 
nately, such studies are often limited either because (1) 
they are purely magnitude limited (resulting in ineffi- 
cient selection of galaxies at the redshifts of interest and 
even fewer galaxies with secure spectroscopic redshifts 
and poorly determined contamination fraction), (2) they 
generally rely on photometric redshifts that are highly 
uncertain at z ~ 2, and/or (3) they are estimated over 
a relatively small number of fields such that cosmic vari- 
ance may be an issue. While purely magnitude limited 
surveys allow one to easily quantify the selection func- 
tion, as we show below, Monte Carlo simulations com- 
bined with accurate spectroscopy can be used to quantify 
even the relatively complicated redshift selection func- 
tions and biases of color-selected samples of high redshift 
galaxies. This "simulation" approach allows one to assess 
a number of systematics (e.g., photometric imprecision, 
perturbation of colors due to line strengths, etc.) and 
their potential effect on the derived LF; these systematic 
effects have been left untreated in previous calculations of 
the L Fs at z ~ 2-3 (jGabasch et al.ll2004lLe Fevre et al.l 
2005), but are nonetheless found t o be important in ac- 
curately computing the LF (e.g.. lAdelberger fc Steidell 
2000; IBouwens et all[200l l27Ml200^ . 

For the past several years, the main focus of our group 
has been to assemble a large sample of galaxies at the 
peak epoch of galaxy formation and black hole growth, 
corresponding to redshifts 1.5 < z < 2.6, in multiple in- 
dependent fields. The selection criteria aim to identify 
actively star-forming galaxies at z ~ 2 with the same 
range in intrinsic UV color a nd extinction as L yman 
break galaxies (LBGs) at z ~ 3 (ISteidel et alJl2003D. The 
color se lection criteria are d escribed in lAdelberger et al.l 
(|2004h : ISteidel et aD (|2004l ). Initial results from the 
survey, including analyses of the star formation rates, 
stellar populations, stellar and dynamical masses, gas- 
phase metallicities, morphologies, outflow properties, 
and clustering are presented in several papers (e.g. , 
Shaplev et alj|2005l; lAdelberger et al . 2005a. b; lErb et al 



logical star formation history (e.g.. ISteidel et al] 
l Adelberger fc Steidell 120001 IBouwens et alj 12004 



2006a|lblld: IReddv et all l2006bL 120051 : IReddv fc Stride! 
2004 ISteidel et al.1 120051 ). With a careful accounting of 
extinction, photometric imprecision, and systematic ef- 
fects caused by observational limitations, rest-frame UV 
selected samples can be used to estimate the cosmo- 
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2006; Bunker et al 1120041 12006^ 

There are primarily two methods by which one can 
attempt to construct "complete" luminosity functions 
that make a reasonable account of all star formation 
at a given epoch. The first method is to observe 
galaxies over as wide a range in wavelengths as possi- 
ble in order to establish a census of all galaxies which 
dominate the star formation ra te density. For exam- 
ple, the union of r est-frame UV (lAdelberger et all 12004 
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(e^ 
1998 



I2004h rest-frame optical (jFranx et al.ll2003 



2004bf). and submillimeter-selected samples 
Hughes et al.1! 1998; Bargc r~et al] 
Chapman et al.ll2005f) should ac- 



Smail et al. 1997 



Blain et al. 2002; 



count for optically-bright galaxies with little to moder- 
ate dust extinction as well as the population of optically- 
faint and heavily reddened galaxies. One can then exam- 
ine the intersection between these various samples and, 
taking into account overlap, compute t he total star for- 
mation rate density (|Reddv et al.ll206~5h . Unfortunately, 
this technique poses several challenging problems, not 
the least of which are the practicality of obtaining multi- 
wavelength data in a large number of uncorrelated fields, 
disparate data quality and photometric depth between 
optical and near-IR images, and the inefficiency of spec- 
troscopically identifying galaxies in near-IR selected sam- 
ples to properly quantify the selection function. 

The second approach, and the one which we adopt in 
this paper, is to estimate sample completeness by way of 
simulations. This method involves simulating many real- 
izations of the intrinsic distribution of galaxy properties 
at high redshift, subjecting these realizations to the same 
photometric methods and selection criteria as applied to 
real data, and adjusting the simulated realizations un- 
til convergence between the expected and observed dis- 
tribution of galaxy properties is achieved. The method 
thus corrects for a large fraction of galaxies that might 
be "missing" from the sample, just as long as some of 
them are spectroscopically observed. The obvious disad- 
vantage of this method is that some (e.g., optically- faint) 
galaxies will never be scattered into our selection window 
and hence we cannot account for such galaxies in our 
analysis. However, multi-wavelength data in several of 
our fields enable us to quantify the magnitude of, and cor- 
rect for, the incompleteness resulting from objects that 
never scatter into our sample. Simulations such as the 
kind presented in this paper become even more important 
at higher redshift (z > 4) where no corresponding multi- 
wavelength data exist to assess the fraction of galaxies 
that are not recovered by color selection. In our case, 
applying the Monte Carlo method to joint photometric 
and spectroscopic samples of high redshift galaxies al- 
lows one to assess the systematic effects of photometric 
scattering and the intrinsic variation in colors due to line 
emission and absorption with unprecedented accuracy. 
In this paper, we take advantage of a large sample of 
spectroscopically-confirmed star-forming galaxies to ar- 
rive at the first completeness-corrected spectroscopic es- 
timate of the UV LF and star formation rate density 
(SFRD) at z ~ 2, computed across the many indepen- 
dent fields of our survey. We extend our results by using 
spectroscopy of Lyman Break galaxies in many new in- 
dependent fields to recompute the UV LF and SFRD at 
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z~ 3. 

While considerable progress in quantifying the cosmic 
star formation history can be achieved by UV obser- 
vations alone, the most robust determination can only 
come from an analysis at multiple wavelengths, where 
systematic effects (e.g., extinction) can be corrected for. 
In addition, assessing the star formation history consis- 
tently at several different wavelengths allows for a useful 
cross-check between results and may reveal any underly- 
ing trends between the star-forming properties of galaxies 
and redshift. In this paper, we combine extensive multi- 
wavelength data in our fields with our spectroscopically- 
derived completeness corrections to measure the rest- 
frame UV, Ha, and infrared luminosity functions at red- 
shifts z ~ 2 — 3. The primary goal of this paper is to 
then use these luminosity functions to evaluate the cos- 
mic star formation history in a consistent manner across 
4 decades of wavelength. 

The outline of this paper is as follows. In § [21 we 
describe the fields of our survey and the color criteria 
used to selected candidate galaxies at z ~ 2. We then 
proceed with a description of the spectroscopic followup 
and quantify the fraction of contaminants, including low 
redshift (z < 1) star forming galaxies and low and high 
redshift AGN and QSOs, within the sample. We con- 
clude §[21 by demonstrating that the redshift distribution 
for the spectroscopic sample is not significantly biased 
when compared with the redshift distribution of all pho- 
tometric candidates at z ~ 2. In §[3j we detail the Monte 
Carlo method used to assess both photometric bias and 
error, the effect of Lya line perturbations on the observed 
rest-UV colors of galaxies, and the procedure used to cor- 
rect our sample for completeness. Our results pertaining 
to the intrinsic Lya equivalent width and reddening dis- 
tributions of 1.9 < z < 3.4 galaxies are discussed in §01 
Results on the rest-frame UV, IR, and Ha LFs are pre- 
sented respectively in § [5j [6l and [71 Lastly, we discuss the 
implications of our results for the luminosity and global 
star formation rate densities in § [HI A flat ACDM cosmol- 
ogy is assumed with Hq = 70 km s _1 Mpc -1 , = 0.7, 
and f2 m = 0.3. 

2. DATA: SAMPLE SELECTION AND SPECTROSCOPY 

2.1. Fields 

Our z ~ 2 survey is being conducted primarily in 
fields chosen for having V < 17.5 mag QSOs with 
redshifts 2.5 < z < 2.8, ideally placed to study the 
correlation between z ~ 2 galaxies and HI and high- 
met allicity (e.g., CIV) abso rbing systems in the IGM 
(see lAdelberger et~a l. 2005b). We have e xtended our sur- 
vey to include the GOODS- North field foickinso n et al.l 
[2001 IGiavalisco et al.ll2004h. encompas s ing th e original 
HDF-North field (Williams jEaI| [1991, HOOQl), and the 
Westphal field (currently encompassed by the large Ex- 
tended Groth Strip survey) to take advantage of the 
multi-wavelength data amassed for these fields. Imag- 
ing was conducted under sim ilar conditions as the z ~ 3 
fields of ISteidel et all (I2003T ). The 14 fields of the z ~ 2 
survey are summarized in Table [H (instrum ents used and 
dates of observation are shown in Table 1 of lSteidel et al.l 
l200l . 

We have expanded significantly the number of spec- 
troscopically confirmed LBGs beyond the or i ginal sam- 
ple from 17 fields presented in ISteidel et al.l (|2003f) . by 



including those LBGs selected in the newer fields of the 
z ~ 2 survey. We have used these new spectroscopic 
redshifts, in additi on to those previously published in 
ISteidel et all (|2003f ) , to re-evaluate the UV LF and SFRD 
at z ~ 3. Fields where we have carried out LBG selection 
are also listed in Table [U 

One of the unique advantages of our analysis is that 
we use a large number of uncorrelated fields (14 and 29 
for the z ~ 2 and z ~ 3 surveys, respectively), combined 
with a large sample of spectroscopic redshifts between 
1.9 < z < 3.4, in order to compute the LF, negating the 
need for uncertain normalization corrections to account 
for clustering and cosmic variance. For example, we find 
evidence for significant large scale s tructure within sev - 
eral fields of the z ~ 2 survey (e.g., ISteidel et all feOOS), 
generally characterized by over-densities in redshift space 
above what would be expected given our redshift selec- 
tion function. By averaging results over many fields well 
distributed throughout the sky, we can estimate the LF 
insensitive to variations in large scale structure, and fur- 
thermore estimate the magnitude of the effect of cosmic 
variance on the results. The total area of all of the inde- 
pendent fields of the z ~ 3 survey is ~ 3200 square ar- 
cmin, or close to a full square degree. The z ~ 2 survey 
area is ~ 1900 square arcmin. Despite the smaller area 
covered by the z ~ 2 survey, there are roughly twice as 
many BX candidates as LBGs given the larger surface 
density of the former. 

2.2. Photometry 

Photo metry was per formed using a modified version of 
FOCAS (lValdeslll982f l. Object detection was done at 1Z 
band, and G — 7Z and U n — G colors were computed by 
applying the 7^. -band isophotal aperture s to the G and 
U n images (see ISteidel et al.ll2003l 120041 for further de- 
tails). The optical images have typical depth of 1Z ~ 27.5 
as measured through a ~ 3" diameter aperture (3 a). 
Field-to-field variations in photometry are dominated by 
systematics due to the different instruments, filter sets, 
and slightly varying observing conditions when the fields 
were imaged. These field-to-field systematics are negli- 
gible compared to measurement errors. We have incor- 
porated some of these effects (e.g., seeing, airmass of the 
observation, CCD response, and filter shape) in comput- 
ing the expected colors of galaxies with known intrinsic 
properties. Modeling all the field-to-field variations in 
photometry rapidly becomes a very complex problem, 
infeasible to resolve within a reasonable time frame. The 
remaining biases (e.g., errors in the zeropoints used) are 
discussed in § 13.21 

2.3. Color Selection 

Even with a priori knowledge of the intrinsic proper- 
ties of all z ~ 2 galaxies, constructing a practical set 
of selection criteria to select all galaxies in any desired 
redshift range and reject all others is an intractable prob- 
lem. One extreme is to select all objects down to a 
given magnitude limit, such as in flux- limited surveys 
of high redshift galaxies, but unfortunately such stud- 
ies suffer from significant amounts of foreground con- 
tamination. Color-selected samples have the advantage 
of allowing one to specifically target a desired redshift 
range while minimizing the number of interlopers. Per- 
haps the most successful of the various color criteria that 
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TABLE 1 
Survey Fields 
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a Right ascension in hours, minutes, and seconds. 
b Declination in degrees, arcminutes, and arcseconds. 
c Number of BX candidates to 1Z = 25.5. 

d Number of spectroscopically confirmed BX candidates with redshifts 1.9 < z < 2.7, excluding those whose spectra 
indicate an AGN/QSO. Note that the total numbers of galaxies, excluding AGN and QSOs, with spectroscopic redshifts 
1.9 < z < 2.7 in the survey (total in all fields is 1288) are larger than the numbers given here since a significant fraction of 
LBG and BM candidates lie at these redshifts. 
c Number of LBG candidates to 1Z = 25.5. 

f Number of spectroscopically confirmed LBG candidates with redshifts 2.7 < z < 3.4, excluding those whose spectra 
indicate an AGN/QSO. Note that the total numbers of galaxies, excluding AGN and QSOs, with spectroscopic redshifts 
2.7 < z < 3.4 in the survey (total in all fields is 1058) are larger than the numbers given here since a small fraction of BX 
candidates lie at these redshifts. 

g Number includes 180 galaxies with 25.5 < 1Z < 26.0. 
h Number includes 10 galaxies with 25.5 < 1Z < 26.0. 

have been designed to select high redshift galaxies is rest- 
frame UV c olor selection, initia lly used to target galax- 
ies at z ~ 3 (Steid eTet al.lll99lf), and extended to higher 
redshifts (e.g., Bouwcns et al. 2005, 2004; B unker et al.l 
I2004I iDickinson et~aTll2004 lYan et al.ll2003h . The suc- 
cess of this technique is partly due to its simplicity in 
that only a few broadband filters are required to assemble 
such samples. Further, at lower redshifts (1.4 < z < 3.5), 
galaxies can be spectroscopically observed and precise 
redshifts can be obtained in a short amount of observing 
time on 8 — 10m class telescopes. Color-selected high red- 
shift galaxy surveys will, as a consequence, have rather 
complex selection functions. The approach described in 
§ [3] allows one to quantify such selection functions with 
relative ease. 

The criteria used to select galaxies with redshifts 1.9 < 
z < 2.7 based on their rest-frame UV colors were de- 
signed to recover objects with intrinsic properties simi- 
lar to those of z ~ 3 Lyman Break Galaxies. The colors 
at z ~ 2 were estimated from spectral synthesis anal- 



ysis of 70 LBGs with broadban d U n GlZJK x photome- 
try and spectrosco pic redshifts (|Adelberger et al.l l2004t 
iSteidel et alj [20041 . Initial spectroscopy of z ~ 2 can- 
didates led to a refinement of the criteria used to select 
galaxies at redshifts 1.9 < z < 2.7 to their present form: 

G-TZ>-0.2 

U n -G>G-TZ + 0.2 

G-TZ<0.2(Un-G)+0A 

U n ~G<G-K+ 1.0, (1) 

termed as "BX" selection (jAdelberger et al.1 12004 
ISteidel et all l2004h. with fluxes in units of AB magni- 



tudes ( Oke fc Gunn 1983). Candidates were selected to 



1Z = 25.5 to ensure a sample of galaxies amenable to 
spectroscopic followup. This limit corresponds to an ab- 
solute magnitude at observed 7?.-band that is 0.6 mag 
fainter at z = 2.2 (the mean redshift of BX candidates 
with z > 1) than at z ~ 3. Additionally, we exclude all 
sources with 72. < 19 that are saturated in our images, 
all of which are stars. The above criteria yielded 10013 
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candidates in the 14 fields, with an average surface den- 
sity of ~ 5 arcmin , uncorrected for contamination from 
objects with redshifts z < 1.4 (see § I2.5[) . The number 
of candidates in each field are summarized in Table [TJ 
Note that we did not select BX gal axies in 15 fields of 
the z ~ 3 survey (jSteidel et al.ll2003T ) since spectroscopy 
of these fields was carried out before the z ~ 2 survey 
began. 

The color criteria used to select LBGs a t red shifts 
2-7 < z < 3.4 are published in ISteidel etlll (|2003t ) and 
are summarized here for convenience: 

G-ft< 1.2 

C/„-G>G- ft +1.0. (2) 

These criteria form the superset of the individual sets of 
criteria for "C", "D ", "M", and " MP" candidate types 
given in Table 4 of ISteidel et~ai] (|2003l) . Hereafter, we 
will refer to all these different candidate types as LBGs. 
Candidates were selected to ft = 25.5, except in the field 
Q1422 where the photometric depth allowed selection of 
candidates to ft — 26.0. The number of z ~ 3 candidates 
in each field are also summarized in Table [TJ Given the 
constraints of the color criteria and the ft — 25.5 spec- 
troscopic limit, the combined BX and LBG samples con- 
stitute ~ 25% of the total 1Z and -Kg-band counts to 
1Z = 25.5 and K S (AB) = 24.4, respectively. 

2.4. Spectroscopic Followup 

The spectros copic followup of can didate s is discussed 
extens ively in ISteidel et al.1 (|2003h and ISteidel et al.1 
(|2004D . Of the 10013 BX candidates, we have targeted 
24% (2382 out of 10013) with spectroscopy, yielding 1711 
rcdshift identifications, or a 72% success rate averaged 
over all fields. As discussed in § I2.6i the spectroscopic 
success rate is primarily determined by the observing 
conditions, and subsequent spectroscopy of spectroscopic 
failures indicates they have similar redshift distribution 
as successes. Similarly, of the 5452 LBG candidates, 
1903 were targeted with 1492 successful redshifts. Fig- 
ure [JJ shows arbitrarily normalized redshift distributions 
for the BX and LBG samples. The mean spectroscopic 
redshifts for the BX and LBG samples, when restricted 
to those objects with z > 1, are (z) = 2.20 ± 0.32 and 
(z) = 2.96 ± 0.26, respectively. Preliminary versions of 
these histograms, along with s ample spectra ofBX galax- 
ies and LBGs, are presented in lSteidel et al.1 (|2003L 120041 ) . 
Table [5] lists the spectroscopic fractions relevant for the 
BX and LBG samples. 

2.5. Interloper Contribution and AGN 

The region of color space defined by BX selection (e.g., 
Figure H|) is also expected to include galaxies outside of 
the targeted redshift range, including st ar forming galax- 
ies at z < 0.2 and stars (see Figure 10 of lAdelberger et aT] 
l2004t ). Spectroscopy shows that there is indeed a subset 
of BX candidates that are interlopers — candidates with 
redshifts z < 1.4 — with a much higher contamination 
rate among candidates with 1Z < 23.5, as indicated in 
Table [3l One can impose a rough magnitude cutoff to 
only consider those candidates with 1Z > 23.5, but this 
would preclude the analysis of the bright-end of the BX 
and LBG luminosity distributions, as well as more de- 
tailed studies of the UV spectra of optically-bright ob- 
jects. Other options to reduce the contamination fraction 




z 



Fig. 1. — Arbitrarily normalized spectroscopic redshift distribu- 
tions of galaxies with z > 1.4 in the BX and LBG samples. The 
total number of galaxies represented here is 2569. 

include using the TZ—K color where the associated bands 
no longer bracket strong spectral breaks for low redshift 
sources . For example, the BzK criteria of iPaddi et al.1 
(|2004bh can be used to reduce the foreground contami- 
nation fraction in color-selected samples. 

The interloper fractions are apt to decrease as the sur- 
vey progresses and we become more adept at excluding 
them from masks based on other multi-wavelength data, 
such as their TZ — K colors. However, until now, we 
have not used any of the techniques discussed above to 
actively discriminate against placing possible interlopers 
on slitmasks; doing so would complicate our ability to 
apply the observed contamination fractions to determine 
the interloper rate among all BX/LBG sources. There- 
fore, the fractions in columns (4) and (7) of Table [3] are 
assumed to represent the overall fraction of interlopers as 
a function of 1Z for the photometric samples. For the BX 
sample, most of the contamination at bright magnitudes 
arises from foreground galaxies. For the LBG sample, 
most of the contamination arises from stars. Applying a 
bright magnitude limit of 1Z = 23.5 will reduce the over- 
all contamination fractions of the BX and LBG samples 
to 9% and 3%, respectively. 

The BX sample also includes a small number of broad- 
lined QSOs and broad and narrow line (a < 2000 km s _1 ) 
AGN whose rest-UV colors are similar to those of high 
redshift star forming galaxies, but which show prominent 
(and in some cases broad) emission lines such as Lya, 
CIV, and NV. The detection rate of such sources is ~ 
2.8% (simi lar to the rate foun d among UV-selected z ~ 
3 galaxies; ISteidel et al.ll2002T ). but is a strong function 
of apparent magnitude where all but two of the objects 
with 7Z < 22.0 and z > 1 are QSOs. The fractions of 
spectroscopically confirmed BXs and LBGs that show 
high ionization UV lines indicative of an AGN or QSO ar e 
listed in Table [2] As discussed in iReddv et all (|2006af l. 
we have found the presence of additional AGN in the 
sample based on either X-ray or IR data. For the analysis 
presented here, we exclude AGN from the sample based 
on the presence of high ionization UV emission lines. The 
effect of including AGN/QSOs in the rest-frame UV LF 
is discussed further in § 18.21 
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TABLE 2 

Spectroscopic and AGN/QSO Fractions of the BX and LBG Samples 





n 




N spec h 


BX 

/spec C 


/ AGN d 


fAPNfz > 1.4) C 




jV apec s 


LBG 

/spec h 


,/agn 


fAOTjtz > 1.4V 


19.0 - 


22.0 


620 


74 


0.12 


0.12 


0.78 


142 


30 


0.21 


0.33 


1.00 


22.0 - 


22.5 


162 


31 


0.19 


0.00 


0.00 


34 


7 


0.21 


0.17 


1.00 


22.5 - 


23.0 


252 


77 


0.31 


0.05 


0.20 


71 


25 


0.35 


0.23 


0.67 


23.0 - 


23.5 


466 


178 


0.38 


0.01 


0.02 


137 


62 


0.45 


0.07 


0.10 


23.5 - 


24.0 


1053 


330 


0.31 


0.03 


0.03 


392 


177 


0.45 


0.02 


0.03 


24.0 - 


24.5 


1894 


511 


0.27 


0.01 


0.01 


881 


398 


0.45 


0.02 


0.03 


24.5 - 


25.0 


2741 


341 


0.12 


< 0.01 


< 0.01 


1617 


442 


0.27 


0.01 


0.01 


25.0 - 


25.5 


2819 


169 


0.06 


0.02 


0.02 


1994 


336 


0.17 


< 0.01 


< 0.01 


25.5 - 


26.0 












180 k 


15 k 


0.08 k 


0.00 k 


0.00 k 


Total 




10007 


1711 


0.17 


0.02 


0.02 


5448 


1492 


0.27 


0.03 


0.03 



a Number of BX candidates. 

b Number of BX candidates with spectroscopic redshifts. 

c Fraction of BX candidates with spectroscopic redshifts. 

d Fraction of AGN/QSOs in BX sample with spectroscopic redshifts. 

e Fraction of AGN/QSO in BX sample with 2; spcc > 1.4. 

f Number of LBG candidates. 

g Number of LBG candidates with spectroscopic redshifts. 

h Fraction of LBG candidates with spectroscopic redshifts. 

1 Fraction of AGN/QSOs in LBG sample with spectroscopic redshifts. 

j Fraction of AGN/QSO in LBG sample with z spcc > 1.4. 

k Numbers are for Q1422 field. 



TABLE 3 

Interloper (z < 1.4) Statistics of the BX and LBG Samples 







BX 






LBG 




n 


AT z > a 


No<z<1.4 b 


/o<z<1.4° 


iVz>o a 


No<z<1.4 b 


/o<z<1.4 C 


19.0 - 22.0 


74 


65 


0.88 


30 


20 


0.67 


22.0 - 22.5 


31 


28 


0.90 


7 


5 


0.71 


22.5 - 23.0 


77 


56 


0.73 


25 


17 


0.68 


23.0 - 23.5 


178 


65 


0.37 


62 


19 


0.31 


23.5 - 24.0 


330 


58 


0.18 


177 


19 


0.11 


24.0 - 24.5 


511 


37 


0.07 


398 


13 


0.03 


24.5 - 25.0 


341 


19 


0.06 


442 


1 


< 0.01 


25.0 - 25.5 


169 


4 


0.02 


336 


2 


< 0.01 


25.5 - 26.0 








15 





< 0.01 


TOTAL 


1711 


332 


0.19 


1492 


96 


0.06 



a Number of sources with spectroscopic redshifts. 
b Number of sources with z < 1.4. 
c Fraction with z < 1.4. 



2.6. Spectroscopic Completeness 

Assessing photometric and spectroscopic completeness 
is a key ingredient in determining the total completeness 
of our survey. The photometric completeness (i.e., the 
fraction of galaxies at redshifts 1.9 < z < 3.4 that satisfy 
either the BX or LBG color selection) is discussed in § [3j 
Here we focus on the extent to which the redshift dis- 
tribution of the spectroscopic sample reflects that of the 
photometric sample as a whole. There are several obser- 
vations that suggest that the redshift selection functions 
for the spectroscopic samples reflect the overall redshift 
selection functions had we obtained spectroscopic red- 
shifts for every single candidate. First, the success of 
measuring redshifts is primarily a function of the weather 
conditions (e.g., cirrus, seeing) at the time of observa- 
tion, with a 90% success rate in the best conditions. Re- 
peat observations of objects for which we were unable to 
secure a redshift initially indicate that the redshift dis- 



tribution of spectroscopic failures is similar to that of 
spectroscopic successes. In other words, our failure to 
measure a redshift is generally not attributable to the 
redshift being far from what one would expect from the 
color selection criteria. 

Second, Figure [5] demonstrates that optical apparent 
magnitude is independent of redshift for the BX and LBG 
samples, keeping in mind that the 1Z — 25.5 limit is ap- 
plied to the photometric (and hence also spectroscopic) 
sample. This is important because if the redshifts of 
objects were correlated with their optical apparent mag- 
nitude, then we might expect the redshift distribution 
of spectroscopically identified candidates to differ from 
candidates in general given that our mask prioritization 
scheme gives more weight to candidates with magnitudes 
in the range 23.5 < H < 24.5 (§E3). Given these results, 
we proceed under the assumption that spectroscopic se- 
lection does not significantly bias the recovered redshift 
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distribution relative to that of the underlying photomet- 
ric sample. 




Fig. 2. — Apparent magnitude versus redshift for spectroscop- 
ically confirmed BX objects (in blue) and LBGs (in red) in the 
redshift range 1.4 < z < 4.0. AGN/QSOs in the BX and LBG sam- 
ples are denoted by the large open circles. The dashed horizontal 
line indicates the self-imposed TZ = 25.5 limit to the photometric 
(and hence spectroscopic) samples. The few objects shown with 
magnitudes fainter than this limit are from the Q1422 field. 



It is instructive to note that, with respect to the red- 
shift distribution, there are two forms of completeness we 
must concern ourselves with. The first is how well the 
redshift selection function for the spectroscopic sample re- 
flects the underlying selection function for the photomet- 
ric sample. We have just argued that the spectroscopic 
and photometric samples must have similar redshift dis- 
tributions. The second is how well the photometric selec- 
tion function reflects the underlying redshift distribution 
of all star-forming galaxies. As discussed in § 15 . 51 the 
number density of galaxies is roughly constant as a func- 
tion of redshift in the redshift ranges 1.9 < z < 2.7 and 
2.7 < z < 3.4. The modulation of this intrinsic (roughly 
constant) redshift distribution into the Gaussian distri- 
bution of Figure [T] can be modeled with great precision 
by way of Monte Carlo simulations, as we demonstrate 
below. Readers who wish to skip directly to the results 
may proceed to § 2] 

3. METHOD: INCOMPLETENESS CORRECTIONS 

A primary aim of this analysis is to connect the ob- 
served properties of BX galaxies and LBGs to the un- 
derlying population of all star-forming galaxies at red- 
shifts z ~ 2 — 3. To this end, we have constructed a 
plausible population of galaxies with a range of redshifts 
(1.4 < z < 4.0), luminosities, and reddening, and deter- 
mined the fraction of these galaxies that would satisfy 
our color criteria. As is typically done, inverting these 
fractions and applying them to the observed counts al- 
lows one to estimate the underlying distribution of galax- 
ies. In this section we discuss in detail the procedure 
used to reconstruct the intrinsic population of galaxies 
at redshifts 1.9 <z < 3.4. 



3.1. Monte Carlo Simulations 

We employed a Monte Carlo approach both to (1) de- 
termine the transformation between the intrinsic proper- 
ties of a galaxy (e.g., its luminosity, reddening, and red- 
shift) and its observed rest-UV colors and (2) quantify 
the effects of photometric errors in their measured rest- 
UV co l ors, similar to the method us ed in Shaplc v" et al.l 
(pOOlj) . lAdelberger fc Steidel pOOOf ). and Steidel et all 



1999). Template galaxies with intrinsic sizes of 0'.'05 to 



0'.'8 and exponential light profiles were convolved with 
the average PSF (typically 1") of the optical images. 
Variations in the light profile used (e.g., exponential, 
de Vaucouleurs) have a negligible effect on the simula- 
tion results; the intrinsic size of the light emit ting region 
(~ 0" 5-0"8 based on HST ACS observations; lLaw et all 
2007) is almost always smaller than the seeing disk. 

The expected rest-UV colors of a galaxy with a par- 
ticular redshift and reddening a re computed by assum- 
ing a iBruzual fe Chariot! (|1996t ) templat e galaxy with 
consta nt star formation for 1 Gyr and a ICalzetti et al.l 
(2000) extinction law. 8 The BX selection criteria were 
designed to select z ~ 2 galaxies with a range of SEDs 
similar to those found for LBGs at higher redshifts 
(jAdelberger et all [2004). Spectral synthesis modeling 
and external multi- wavelength information indicates that 
most UV-selected z~2 — 3 galaxies can be described by 
long duration (> 100 Myr) starbursts and the constant 
star formation model described above should reproduce 
this behavior to the exte nt required by the simulations 
(e.g- JShaplev et al" , 2005). In particular, the rest-UV col- 
ors of galaxies are essentially constant after 10 s years of 
star formation, onc e the m ix of O and B stars stabilizes 9 . 
The ICalzetti et al.l ()2000D reddening law reproduces the 
average expected star formation rates of z ~ 2 — 3 galax- 
ies based on extinction free stacked X-ray and radio es- 
timates fe.g.. lReddv fc Steidell l2004h and further repro- 
duces the average dust obscuration of galaxies with bolo- 
metric luminosities in the range 10 11 < -^boi S, 10 12 2 Lw 
where the bulk of our sample lies (|Reddv et al J [2006b) . 
The use of a constant star forming model and the Calzetti 
reddening law should therefore adequately parameterize 
the SEDs of most optically-bright star-forming galaxies 
at z ~ 2 — 3. An advantage of spectroscopic followup of 
photometrically selected BX galaxies and LBGs is that 
we can also constrain the effects of IGM opacity and 
Lya absorption/emission (§ I3.3|) . both of which are red- 
shift dependent. All of these perturbing effects will re- 
sult in a wide range of spectral shapes and should ac- 
c ount for any galaxies that are not exactly described by 
a lCalzetti et al.l ()2000f ) attenuated constant star forming 
SED. 

A large distribution of galaxy colors was then com- 
puted assuming a particular luminosity function (LF) 
and the observed E(B — V) distribution for spectroscopi- 
cally confirmed galaxies. Small variations in the assumed 

8 Note that because our selection, and hence simula- 
tions, are concerned with the rest-UV colors, adopting a 



IMarasto n et al .l II2006T1 model (where most of the difference with 
the IBruzual fc Chariot! | Tl996l) model is in the rest-optical) should 
minimally affect our results. 

9 There is considerable leeway in the best-fit star formation 
histories for the optical/IR SEDs of UV-selected z ~ 2 galaxies, 
but external constraints point to bur st timescales of > 100 Myr 
llShaplev et a l. 2005; Erb et al. 2006c). 
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Schechter parameters of the LF do little to change the 
results, since our main goal is to sufficiently populate 
rcdshift space and rest-UV color space with a realistic 
distribution of objects. The results are also insensitive 
to small variations in the assumed E(B — V) distribution 
as long as the range of E(B — V) chosen reflects that 
expected for the galaxies. A by-product of the luminos- 
ity function analysis is that we also compute the best-fit 
underlying E[B — V) distribution. The validity of the 
assumed LF and E{B — V) distributions can be tested 
by comparing with the inferred LF and E(B — V) dis- 
tributions. Significant differences between the assumed 
and inferred distributions imply that the initial assump- 
tions for the LF and E(B — V) distribution were differ- 
ent from their true values. The colors were corrected for 
o pacity due to t he intergalactic medium (IGM) assuming 
a lMadaul (|1995[ ) model, and corrected for filter and CCD 
responses and airmasses appropriate for the each field of 
the survey. 

The intrinsic rest-UV colors are randomly assigned to 
simulated galaxies that are then added to the images in 
increments of 200 galaxies at a time. This ensures that 
the image including all added (simulated) galaxies has 
confusion statistics similar to the observed image, since 
this will affect the photometric uncertainties and sys- 
tematics due to blending. We then attempt to recover 
these simulated galaxies using the same software used to 
recover the real data, and record whether a simulated 
galaxy is detected and what its observed magnitude and 
colors are. We repeated this procedure until approxi- 
mately 2 x 10 5 simulated galaxies were added to each of 
the U n , G, and 1Z images of each field. This large number 
of simulated galaxies is necessary in order to sufficiently 
populate each bin of luminosity, reddening, and redshift. 

The end product of the simulations are sets of trans- 
formations for each field that give the probabilities that 
galaxies with intrinsic luminosities (£'), reddenings {E'), 
and redshifts (2/) will be observed to have luminosities L, 
reddenings E, and redshifts z (or alternatively, the prob- 
abilities that galaxies with true properties L'E'z' will be 
measured with a particular set of rest-UV colors). 

3.2. Photometric Uncertainties 

We have used the results of the Monte Carlo simu- 
lations (§ 13. 1|) to estimate the photometric errors and 
determine optimal bin sizes for subsequent analysis. For 
each simulated galaxy that is detected, we have recorded 
the true and measured rest-UV colors. As the uncertain- 
ties may vary depending on magnitude or color, we have 
binned the detected galaxies in magnitude and color for 
each field and have only considered galaxies that would 
be detected as candidates since these are the only objects 
that are relevant to our analysis. We used bin sizes of 
0.5 mag in 1Z and 0.2 mag in U n — G and G — 1Z color to 
determine the uncertainties in the recovered magnitudes 
and colors of objects in each field. Systematic bias in 
the G — 1Z color was estimated by computing the quan- 
tity A[G — 7Z] = (G — lZ) mcas — (G — 1Z) t ruc which was 
typically < 0.04 mag with uncertainty estimated to be 
<j(A[G — 1Z]) ~ 0.09 mag. The typical random uncertain- 
ties in U n — G and 1Z are ~ 0.15 mag and ~ 0.13 mag, 
respectively. These quantities were determine d usin g 
the same me t hod a s pre sented in IShaplev et all (|2005f l. 
iSteidel et al.l ([2003), and IShaplev et al.l (|200lD The un- 



certainties were ge nerally larger for objects faint in 1Z 
(|Steidel et al.ll2003f l The field-to-field results were con- 
sistent with each other (i.e., the typical biases and un- 
certainties from field-to- field were within 0.1 mag of each 
other). The photometric errors are slightly smaller in size 
than the bin sizes (0.2 mag) used to estimate the redden- 
ing and luminosity distribution. A more refined method 
discussed in § 13.4.21 will correct for any systematic scat- 
tering of objects into adjacent bins due to photometric 
error and/or Lya perturbations to the colors. 

3.3. Lya Equivalent Width (Wi,y a ) Distribution 

The presence of Lyman alpha absorption and/or emis- 
sion can perturb the observed rest-UV colors of z ~ 2 — 3 
galaxies by up to 0.75 mag depending on the redshift and 
intrinsic (rest-frame) Lya equivalent width W-Lya- To 
investigate these effects, we measured the W^ya for 414 
spectroscopically confirmed BX galaxies with redshifts 
1.9 < z < 2.7 in 7 different fields of the BX survey. The 
resulting distribution for all galaxies with 1.9 < z < 2.7 
is shown in Figure [3K, with the characteristic asymmet- 
ric shape. The photometric scattering probability asso- 
ciated with this Wi,y a distribution is shown in Figure 2J 
The two shaded "zones" in Figure H] reflect the redshift 
ranges where Lya falls within the U n and G-bands. This 
figure demonstrates how galaxies that are targeted by 
the BX criteria can be shifted out of the BX selection 
window due to Lya emission or absorption. 
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Fig. 3. — (a) Rest frame Lya equivalent width (W^ ya ) distri- 
bution for 482 spectroscopically observed z ~ 2 galaxies. Panels 
(b)-(d) show the Wi^ ya distribution for subsets in redshift. We 
use the convention that W-^ya > implies emission. The distri- 
butions are absorption dominated in all cases. The £ values in- 
dicate the probability that the distributions are drawn from the 
same parent population as the Wx,y a distribution for galaxies at 
2.17 < z < 2.48, the redshift range where Lya does not affect the 
rest-UV colors. 

Since our ultimate goal is to determine how the WL yct 
distribution perturbs the intrinsic colors of galaxies (i.e., 
the colors we would measure in the absence of absorp- 
tion and/or emission), we must first determine whether 
the measured Wh ya distribution reflects the intrinsic 
distribution for the parent population of galaxies. In 
other words, we must check if our color selection crite- 
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Fig. 4. — Perturbation of U n G1Z colors from Ly« absorption and 
emission. The trapezoid is the BX selection window defined by 
Equation [T] The U n GTZ colors of a template galaxy with constant 
star formation for > 100 Myr (after which the UV colors are es- 
sentially constant) and E(B — V) = 0.13 (the mean for the z ~ 2 
sample) assuming a Calzctti ct al. (2000) extinction law is shown 
by the solid curve, proceeding from redshift z = 1 to 3. The lower 
and upper shaded regions correspond to redshift ranges where the 
Lya line falls in the Un and G-bands, respectively. In the absence 
of photometric errors and assuming all galaxies can be described 
by the SED assumed here, galaxies with redshifts 1.68 < z < 2.17 
and 2.48 < z < 2.93 will fall in the dark gray regions with a prob- 
ability of 64% based on the H^Lya distribution in Figure [3] The 
medium and light gray regions correspond to scattering probabil- 
ities of 30% and 6%, respectively. Arrows labeled "abs" indicate 
the direction in which the colors will be perturbed with increasing 
Lya absorption. 

ria introduces significant biases into the measured W^ ya 
distribution. We can begin by examining some char- 
acteristics of the measured Wh ya for BX galaxies and 
LBGs, summarized in Table |4] The BX distribution has 
(Why a ) ~ — 1 A, somewhat lower than the mean Why a 
for LBGs. While the measurements of Wh ya for individ- 
ual galaxies may be highly uncertain, the difference in 
the average Wh ya suggests that the high redshift (LBG) 
population has a higher incidence of Lya in emission than 
the low redshift (BX) population. This disparity between 
the lower and higher redshift populations can be better 
appreciated by examining column (3) of Table H that 
shows that the fraction of galaxies with Wh ya > 20 A 
(/20) is almost twice as high among LBGs as it is for 
BX galaxies. 

The change in /20 is even more apparent when we 
consider BX galaxies in different redshift ranges: /20 for 
BX galaxies with redshifts between z « 2.5 and z = 2.7 
is twice that of BX galaxies with redshifts between z = 
1.9 and z « 2.5 (Tabled}. These results suggest that 
star-forming galaxies exhibit a higher incidence of Lya 
emission at higher redshifts, but is this trend intrinsic to 
high redshift galaxies, or is it introduced as a result of 
color selection bias, as Figure [4] suggests? 

We can test for systematics induced by the color cri- 
teria by examining /20 for BX galaxies at redshifts 
2.17 < z < 2.48, where Lya lies outside the U and 
G-bands. These galaxies have a similar /20 to that of 
z < 2.17 galaxies (Table SJ, implying that the frac- 
tion of absorption versus emission line systems culled 
by the BX criteria is similar between the z < 2.17 and 



TABLE 4 
Measured WL ya Distributions 



Sample 


(W Lya )* 


f(W hya ) > 20 A b 


BX (ALL: 1.90 < z < 2.70) 


-1 A 


0.12 


BX (1.90 < z < 2.17) 


-1 A 


0.08 


BX (2.17 < z < 2.48) 


-2 A 


0.11 


BX (2.48 < z < 2.70) 


2 A 


0.20 


LBG (2.70 < z < 3.4) 


9 A 


0.23 



a Mean rest-frame Wi^ ya . 
b Fraction with W Lya > 20. 



the 2.17 < z < 2.48 samples (this assumes that there 
is little evolution in the Wh ya distribution between the 
z < 2.17 and 2.17 < z < 2.48 subsamples). Focusing 
on the high redshift subsample with 2.48 < z < 2.70, 
Figure [4] suggests that these galaxies are more likely to 
satisfy the BX criteria if they have Lya in absorption, 
yet their /20 is similar to that of LBGs. In other words, 
the 2.48 < z < 2.70 subsample has an /20 value that 
does not indicate a preferential selection of absorption 
over emission line galaxies relative to that of the lower 
redshift subsamples. Rather, the /20 value is larger than 
those for the lower redshift subsamples and is similar to 
that of the LBGs. These conclusions are supported by a 
Kolmogorov-Smirnov test, the results of which are sum- 
marized in Figure [3] Namely, £ in the figure indicates 
the probability that the Wh ya distributions for the to- 
tal sample, the sample with 1.90 < z < 2.17, and the 
sample with 2.48 < z < 2.70, are drawn from the same 
parent population as the sample with 2.17 < z < 2.48, 
where Lya does not effect the U n GTZ colors. Galaxies 
with 2.48 < z < 2.70 have a Wh ya distribution that de- 
viates significantly from the one at 2.17 < z < 2.48. 

Assuming that the UV properties of galaxies are inde- 
pendent of their Lya line profiles would then suggest that 
the BX color criteria do not significantly modulate the 
intrinsic Wh ya distribution of z ~ 2 galaxies. 10 For the 
purposes of our simulations, we make the approximation 
that the observed Wh ya distribution for BX galaxies can 
be applied to our simulated galaxies to obtain the aver- 
age perturbation of their rest-UV colors. 

Since the Lya line falls in the G-band for galaxies in 
the entire redshift range 2.7 < z < 3.4, we cannot ex- 
amine trends in the Wh ya distribution for the LBGs in 
the same way we did for the BXs. However, in § 14.11 
we justify why the Wh ya distribution for LBGs should 
approximately reflect the intrinsic Wh ya distribution for 
z ~ 3 galaxies. Also in § 14.11 the incompleteness correc- 
tions are used to test whether our initial assumptions of 
the Wh y a distributions are correct. 

3.4. Quantifying Incompleteness 
3.4.1. Effective Volume (V e $) Method 

The fraction of galaxies with a given set of binned prop- 
erties that satisfy the color criteria can be computed di- 
rectly from the results of the Monte Carlo simulations. 

10 IShaplev et al.l l|2003D have demonstrated that WLya is in fact 
dependent upon the rest-frame UV colors and magnitudes of galax- 
ies. However, the small biases that these trends may have on the 
observed Wi, ya do not have a significant impact on the derived LFs 
at z ~ 2 and z ~ 3. 
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These binned properties might be the optical luminosity 
(L), redshift (z), and reddening (E(B — V)) of a galaxy. 
Under the assumption that these properties are indepen- 
dent of each other, and if we let the indices i, j, and 
k run over the range of values of L, z, and E(B — V), 
then the true number of galaxies in the ijkth bin can be 
approximated as 



Bins of Intrinsic Properties 



true 
l ijk 



l ijk 



/Pijk 



(3) 



where pijk are the mean probabilities that a galaxy in 
the ijkth bi n is (a) detected and (b) satisfies the color 
criteria (e.g., lAdelbergerir2002t ). These probabilities pijk 
are simply 



Pijk 



1 ™ 
— ^Pvt 



(4) 



where Pijkn is the probability that the nth simulated 
galaxy in the ijkth bin will be detected as a candidate, 
and riijk is the total number of simulated galaxies in the 
ijkth bin. The values Pijkn take into account the prob- 
ability that the colors of the nth simulated galaxy will 
be perturbed by the W^ya distribution of Figure [3] and 
still be selected as a BX object. They also fold in the 
probability that a non-candidate simulated galaxy will 
fall in the BX selection window. Dividing by riijk nor- 
malizes the mean probabilities pijk and accounts for both 
the fraction of galaxies whose photometric errors scatter 
them out of the BX selection window and galaxies that 
are not detected in the simulations. If the true comov- 
ing volume corresponding to the j th bin in redshift is Vj , 
then the effective volume associated with the jth bin in 



V, 



eff 



ik 



Vi > Pijk 



Vj x & 



(5) 



where £j are commonly referred to as "completeness co- 
efficients" : 



(G) 



The photometric properties of each field are unique due 
to differences in the observing conditions, and this will 
affect the computed £j . We can then determine the com- 
pleteness coefficients for each field and then perform a 
weighted-average of them (i.e., weighted according to the 
field size) to obtain mean completeness coefficients, £j. 

3.4.2. Maximum Likelihood Method (V\ik) 

While the procedure just described can be used to 
make an initial guess as to the shape of the reddening 
and luminosity distributions, it can lead to spurious re- 
sults, particularly for objects whose true colors are such 
that they lie outside of or close to the edges of the BX 
selection window. Equation [3] is approximately true only 
if the average measured properties of a galaxy are the 
same as the true (simulated) properties, and this will 
certainly not be the case for galaxies that are preferen- 
tially scattered into the BX window due to photometric 
errors or the pres ence of Lya absorption/emission (e.g., 
lAdelbergerl 120021 ). The approach described above will 
also not take into account photometric bias and the pref- 
erential scattering of objects from one bin to another if 



Po'O'C 




Bins of Observed Properties 

Fig. 5. — Cartoon illustration of how the probability that a 
galaxy with intrinsic (true) properties L 1 E' z 1 may not have a one- 
to-one correspondence with bins of observed (measured) properties 
LEz (or measured colors ugr). 

the bin sizes are c omparable to (or smaller than) the 
photometric errors (|Adelbergerll2002f ). 

Figure [5] further illustrates these issues. In the sim- 
plest case, galaxies that fall within a particular bin of 
true properties (say, L' E' z' ) will, on average, have mea- 
sured properties corresponding to bin LqEqZq. In this 
case, we can use the approach of § 13.4.11 (i.e., the V b q 
method) to simply divide the observed number of galax- 
ies in bin LqEqZq by the probability that galaxies in bin 
LqEqZq will be observed in bin LqEqzo (call that prob- 
ability po'O'O'^ooo), as shown by the leftmost arrow in 
Figure [5j However, we can point to several examples 
that suggest that there may not be a one-to-one corre- 
spondence between bins of intrinsic and observed prop- 
erties, as illustrated by the remaining arrows in Figure[5j 
First, in order to accurately compute the luminosity and 
reddening distributions at z ~ 2, we cannot make our 
bin sizes much larger than the typical photometric errors 
since the observed range of U n — G and G — 1Z colors for 
galaxies at a single redshift (< 0.8 mag) is only < 4 times 
the typical photometric error in U n —G and G—1Z. There- 
fore, galaxies that ought to fall within a particular bin of 
measured properties will be scattered into adjacent bins. 
This would not be a problem if each bin of measured 
properties gained and lost an equal number of galax- 
ies, but since the luminosity and reddening distributions 
are peaked, photometric errors will scatter galaxies away 
from the peak and into the wings of the distributions. 
Second, the distribution of errors in colors is not sym- 
metric with respect to the true values such that there is 
a systematic tendency to sca tter galaxies into red der bins 
more often than bluer ones (jSteidel et al.ll2003h . Third, 
the presence of Lya absorption in a galaxy's spectrum 
will, depending on the redshift, cause us to overestimate 
the reddening. Finally, there will be some galaxies whose 
true properties are such that on average they lie outside 
the selection windows, and only get scattered into the 
sample because of photometric errors, such as might be 
the case for galaxies lying close to the color selection 
boundaries. 

Because of these systematic effects, the number of 
galaxies within a particular bin of measured properties 
will be some weighted combination of the numbers of 
galaxies within intrinsic bins that contribute to that mea- 
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sured bin of properties. The weights are simply the 
"transitional" or scattering probabilities indicated by the 
arrows in Figure and these can be determined from the 
Monte Carlo simulations. Formally, the number of galax- 
ies we expect to observed in a bin of measured luminosity 
L, reddening E, and redshift z within a field of size AQ 
is 



n{L, E, z) 



dV 



fi J dL'dE'dz'f(L')g(E')h{z')p L , E/z ,^ LEz 



dQdz 
dV 



dQdz 



, (7) 



where y, is related to the total comoving number density 
of galaxies; pl'E'z'-»lez is the (transitional) probability 
that galaxies with intrinsic L'E'z' will be measured to 
have LEz; and f(L'), g(E'), and h(z') are the intrinsic 
distributions of luminosity, reddening, and redshift, re- 
spectively, normalized such that 



dLf(L) = 



dEg(E) 
dzh(z) = 



1 



(8) 



( Adelberger 2002). Our goal is to determine the intrinsic 
distributions f{L), g(E), and h(z), but inverting Eq. [7] 
to solve for these distributions is intractable. One al- 
ternative is to compute the likelihood (£) of observing 
our data, which is expressed as a list of galaxies with 
observed LiEiZi, for a given set of fgh distributions: 



C(LiEiZi) oc exp 



-/iAf2 J dLdEdzn(L,E,z 



dV 



dildz _ 
Hn(LiE iZi %9) 



The discrete form of Eq. [51 extended to incorporate each 
of / different fields, can be expressed as 



£(n ijk i) oc exp 



ijkl 



I'ijkl 



n 

ijkl 



'ijkl 



(10) 



where riijki is the mean number of galaxies in the i bin 
of luminosity, j th bin of reddening, and fc th bin of redshift 
in the I th field that the assumed values of fi, gj, and hk 
imply; and na ki is the observed number of galaxies in the 
same bin (e.g.. lAdelbergerj|2002D . The discrete version of 
Eq.0is 



^ijkl 



= fiAQi V fi'9j'hk'Vk>pi,i',jr,k' 



i'j'k' 



*ijk > 



(11) 



where Aili is the size of the I field, Vk> is the comov- 
ing volume in Mpc 3 arcmin -2 corresponding to bin k' in 
redshift, and Pi,v,j' k'—*ijk is the probability that a galaxy 
in the I th field in the i'j'k' bin of luminosity, reddening, 
and redshift, will have measured properties correspond- 
ing to bin ijk. Assuming that the data quality does not 
vary significantly from field-to-field, we can simplify the 
probabilities such that 

Pi'j'k'^ijk = ^ MllPli'j'k'-njk/ ^^Afli. (12) 



Maximizing the likelihood as expressed in Eq.[10]is equiv- 
alent to minimizing 

- m£ oc fjjjk - riijk Innyfe, (13) 

ijk ijk 

and is more amenable to computation than Eq. 1101 

3.4.3. Implementation of the Maximum- Likelihood Method 

We first used the Monte Carlo simulations to determine 
the transitional probabilities that relate the true lumi- 
nosities, reddenings, and redshifts of galaxies to their ob- 
served rest-UV colors. Following the discussion of § I3.3( 
the colors of galaxies were perturbed by randomly as- 
signing a Wiya according to the distributions shown in 
Figure [3] for t he BX sample and th e distribution shown 
in Figure 8 of lShaplev et all (|2003l ) for LBGs (see also 
Figure [5]). We took advantage of both the U n — G and 
G—1Z colors in our analysis of the z ~ 2 sample to provide 
more stringent constraints on the E{B — V) distribution, 
something that was not possible at z ~ 3 where most 
galaxies only had limits in U n either due to severe blan- 
keting by the Lya forest or the suppression of continuum 
flux shortward of the Lyman limit. 

Figure [5] is useful in visualizing the transitional prob- 
abilities, in this case for the BX sample, where we show 
the relative probability distribution for galaxies between 
1.0 < z < 3.0 to be selected by the BX criteria. The 
probability distribution is weighted by the incidence of 
galaxies with intrinsic colors as determined from the LF 
and E{B — V) distributions assumed in computing the 
transitional probabilities. This distribution reflects both 
photometric error and Lya perturbation of the expected 
rest-UV colors. One noticeable feature of Figure [6] is 
the divergent behavior of the selection function for low 
(z < 2.0) and high [z > 2.7) redshift galaxies, where 
higher redshift galaxies have redder U n — G colors for 
a given SED. This can be understood, in part, by ex- 
amining Figure 21 If z ~ 2 galaxies can be reasonably 
described by the SED and reddening assumed above then 
we would expect that galaxies with z > 2.7 would only be 
scattered into the BX window if there were large changes 
in their colors, either due to photometric errors or Lya 
perturbation. First, we find no evidence that photomet- 
ric errors increase for galaxies at higher redshifts. Sec- 
ond, the (l + z) dependence of the observed Wh ya will re- 
sult in a larger color change (for a fixed rest-frame W\^ ya ) 
for higher redshift galaxies than for lower redshift galax- 
ies, such that the scattering probability distribution cov- 
ers a larger area in color space, making it less likely for a 
particular source to fall within the BX selection window. 
Finally, the U n — G color changes more rapidly for higher 
redshift galaxies where Lya forest absorption begins to 
increasingly affect the f/ n -band. All of these effects could 
explain the relatively small number of z > 2.7 galaxies 
singled out with the BX criteria. The advantage of rest- 
UV selection is that the drop off in BX efficiency for 
z > 2.7 can be compensated for by adopting the z ~ 3 
LBG criteria whose selection function begins to rise for 
z > 2.7 and which use exactly the same filt er set, negat- 
ing t he need for additional observations (jSteidel et alJ 
12001 . 

Unlike the z > 2.7 galaxies discussed above, 1.0 < z < 
2.0 galaxies are crowded into a narrower region of color 
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Fig. 6. — Relative probability distribution for galaxies with in- 
trinsic colors (U n — G)truc and (G — 7?,)truc to be detected and 
selected as BX objects (solid line same as in Figure [4j. The dis- 
tribution is weighted according to the incidence of galaxies with 
a particular set of intrinsic colors as determined from the LF and 
E(B — V) distributions used to compute the transformation be- 
tween intrinsic and observed colors. The distribution is non-zero 
exterior to the BX window (trapezoid) as a result of photometric 
error and Lyo line perturbations of the colors. Galaxies with ex- 
pected (or intrinsic) Un — G colors bluer than required to s atisf y BX 
criteria are particularly prone to selection as discussed in § 13.31 The 
region between the white curves denotes the swath of color space 
where galaxies with redshifts 2.17 < z < 2.48 are expected to lie. 
These galaxies' colors are unaffected by Lya line perturbations. 

space as is evident from Figured) Small variations in col- 
ors as a result of photometric errors or Lya absorption 
can shift a large numbers of such galaxies into the BX 
selection window. This effect can be viewed in Figure [5] 
where there is a high relative probability for galaxies with 
bl ue U n - G colors (the "BM" galaxies; e.g., Figure 10 
of lAdelberger et aT1l2004h to satisfy BX selection, partly 
due to the effect of Lya absorption in these systems (cf., 
Figure [3)3). The highest density region in this figure (be- 
tween the two white curves of Figure [6]) occurs in the 
same color space expected to be occupied by galaxies at 
redshifts where the Lya line does not affect the U n GlZ 
colors (2.17 < z < 2.48). Figure [6] also demonstrates the 
fallacy of the assumption in Equation [3] where the true 
and observed rest-UV colors may be significantly and, 
more importantly, systematically different for galaxies ly- 
ing in particular regions of color space. Figure[6]is meant 
to be purely illustrative and, in reality (as in our simu- 
lations), the probability distribution will be "smeared" 
out when one considers galaxies with a range of spectral 
shapes. 

The effects of IGM opacity, Lya absorption/emission, 
and photometric error (§ 13 . 113.31 13 . 2[) imply that simple 
boxcar approximations to the redshift selection function 
(even in photometric surveys) are unrealistic, irrespec- 
tive of the wavelengths used to select galaxies. The ad- 
vantage of our combined Monte Carlo, photometric, and 
spectroscopic approach is that even complicated selection 
functions can be quantified relatively easily and thus be 
corrected for in the final analysis. 

In practice, minimizing Eq. [13] while simultaneously 
varying the intrinsic distributions of luminosity, redden- 
ing, and redshift (fgh) can lead to spurious results given 



the large parameter space and possibility of numerous lo- 
cal minima in likelihood space. A reasonable approach is 
to then make some simplifying assumptions, such as fix- 
ing the redshift distribution to be constant and assuming 
an LF computed using the method of § 13.4.11 One can 
then minimize Eq. 1131 with respect to the distribution of 
spectral shapes (g) as parameterized by the E(B — V) 
color excess, using values of Pi'j>k'->ijk relevant for the 
spectroscopic sample. In other words, pi'j>k'^ijk will 
give the probability that a galaxy with true properties in 
the i'j'fc' th bin will be measured with properties in the 
ijk th bin and be spectroscopically observed. The proba- 
bility that a candidate lying within a particular bin of 1Z 
magnitude will be spectroscopically observed is approxi- 
mated using the spectroscopic fractions listed in Table [2] 
These spectroscopic fractions are then multiplied by the 
probability that an object is a star- forming galaxy (i.e., 
not an AGN/QSO) using the AGN/QSO fractions in the 
relevant magnitude range (Table [2]) . At this stage, we 
must rely on the spectroscopic sample since we can only 
estimate E(B — V) for galaxies with redshifts. Then, 
keeping g fixed to the best-fit E(B — V) distribution, we 
take advantage of the full photometric sample to mini- 
mize Eq. Q2] with respect to the luminosity distribution 
(/). The revised estimate of / can then be held fixed to 
refine our estimate of g. The process goes through several 
iterations where at the last stage we vary /, g, and h si- 
multaneously. The results of this procedure indicate that 
our initial assumption of a constant redshift distribution 
(i.e., number of galaxies in each of the redshift bins is 
roughly constant) is a reasonable one to make. The red- 
shift distributions predicted for BX galaxies and LBGs 
given the maximum-likelihood fgh distributions are ex- 
cellent matches to the observed redshift distributions of 
BX galaxies and LBGs, as shown in Figure [7] 



LBG 





Fig. 7. — Expected redshift distributions (lines) given our best- 
fit reddening and luminosity distributions, compared with the ob- 
served redshift distributions of BX galaxies (left panel) and LBGs 
(right panel), indicated by the shaded histograms. 



Uncertainties in the luminosity and E{B — V) distri- 
bution were estimated by generating many fake realiza- 
tions of our observed data from the catalogs of simulated 
galaxies, and recomputing the best-fit fgh. The disper- 
sion in measurements of fgh are taken to be the 1 a 
errors. It is important to note that the errors in our es- 
timates are due to a combination of Poisson noise and 
field-to-field dispersion. Unlike all other previous esti- 
mates of the z ~ 2 — 3 UV LF, our determination incor- 
porates the largest spectroscopic sample of galaxies at 
these redshifts and automatically takes into account the 
systematic effects mentioned in § 13.4.21 
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Fig. 8. — Comparisons between expected and observed tVLya distributions for different assumptions of the intrinsic W\^ ya distributions 
between redshifts 2.7 < z < 3.4 (top panels) and 1.9 < z < 2.7 (bottom panels). The assumed intrinsic distributions are denoted by green 
rectangles in the left and middle panels. The assumed intrinsic distributions are the same as the observed distributions (blue rectangles) 
in the right panels. The expected and observed distributions are indicated in red and blue, respectively, in all six panels. The length of the 
bars represent the dispersion in values of the tVLya distribution derived assuming many realizations of the LF and E(B — V) distribution 
(see text). The parameter £ denotes the likelihood that the observed and expected distributions are drawn from the same distribution. We 
use the convention that WLya > implies Lyct emission. 



4. RESULTS: INTRINSIC W LYa AND E(B - V) 
DISTRIBUTIONS 

4.1. Validity of Assumed Whya Distributions 

An important question is whether the distribution of 
Lya emission and absorption profiles of galaxies changes 
as a function of redshift. Such trends with redshift may 
indicate fundamental differences in the ISM of galax- 
ies and/or changing large-scale environments as a func- 
tion of redshift. Can we do better job of determin- 
ing whether the intrinsic W^ya distribution of galaxies 
changes as a function of redshift? Ideally, we would have 
liked to include the W^ ya distribution as another free 
parameter in the maximum-likelihood method discussed 
in § I3.4.2|3.4~3l so that instead of maximizing three func- 
tions (fgh), we would be maximizing four. However, this 
would needlessly complicate our ability to determine the 
maximum-likelihood fgh distributions, especially since 



the luminosity distribution (/) is insensitive to small 
changes in the Wi^ya distribution. As a compromise, we 
can investigate how different assumptions of the intrin- 
sic W-Lya distributions of galaxies affect the distributions 
that we expect to measure. 

Figured] illustrates how the color selection criteria can 
modulate the observed WLya distribution of galaxies, 
such that the observed distribution may be different than 
the intrinsic distribution. The Monte Carlo simulations 
discussed in § [3] allow us to directly compare the mea- 
sured (observed) W-Lya distributions for the BX and LBG 
samples with those expected based on the transitional 
probabilities. The results of this comparison are sum- 
marized in Figure [8l which shows W^ ya for various as- 
sumptions of the input W\^ ya distribution. We consider 
three cases. In the first case, we assume that the intrinsic 
W^Lya distribution of galaxies at z ~ 3 is identical to that 
measured at z ~ 2. In the second case, we assume that 
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the intrinsic W^ ya distribution at z ~ 3 is an equally 
weighted combinations of the measured distributions at 
z ~ 2 and z ~ 3. In the third case, we assume that the 
intrinsic WL yQ distribution at z ~ 3 is identical to the 
measured distribution at z ~ 3. Three analogous cases 
are considered for the z ~ 2 sample. 

For example, the top left panel of Figure [5] shows the 
W-Lya distribution for 2.7 < z < 3.4 galaxies that one 
would expect (red rectangles) if the intrinsic distribu- 
tion at 2.7 < z < 3.4 is identical to the measured W^ ya 
distribution for lower redshift (1.9 < z < 2.7) galaxies 
(green rectangles, labeled "Obs z ~ 2"). The validity of 
assuming a particular intrinsic distribution can be tested 
by comparing the expected distribution (red rectangles) 
with the actual measured distribution (blue rectangles). 
The vertical sizes of the rectangles for the observed dis- 
tributions reflect Poisson errors. Uncertainties in the ex- 
pected distributions (red rectangles) are computed by 
constructing many samples of galaxies drawn randomly 
from the maximum-likelihood luminosity and E(B — V) 
distributions (§ I4.21 § IS")), and fixing the vertical bar size 
to the dispersion in the Wi^ ya distributions measured for 
each of these simulated samples. 

The top left panel of Figure [5] shows that assuming 
an intrinsic distribution of WL ya for galaxies at redshifts 
2.7 < z < 3.4 that is identical to the measured distri- 
bution of W^ya for BX-selected (1.9 < z < 2.7) galaxies 
results in an expected distribution at 2.7 < z < 3.4 that 
deviates significantly from the distribution that we ac- 
tually measured. In this case, the expected distribution 
exhibits a larger fraction of galaxies with absorption and 
a deficit of emission-line galaxies when compared with 
the measured W^ ya distribution. Therefore, the intrin- 
sic WLya f° r 2.7 < z < 3.4 galaxies must have lower and 
higher fractions, respectively, of absorption and emission- 
line galaxies than what is observed among lower redshift 
(1.9 < z < 2.7) galaxies. The bottom left panel of Fig- 
ure [5] tells a similar story. Assuming 1.9 < z < 2.7 galax- 
ies have an intrinsic WLya identical to that measured for 
LBGs (2.7 < z < 3.4) results in an expected distribution 
for 1.9 < z < 2.7 galaxies that has a lower frequency of 
absorption-line systems than what is actually observed 
at 1.9 < z < 2.7. Therefore, the intrinsic WLya distri- 
bution for 1.9 < z < 2.7 galaxies must include a larger 
fraction of galaxies with Lya in absorption than what is 
observed for higher redshift (2.7 < z < 3.4) galaxies. 

The middle panels of Figure [8] show what happens 
if we assume that the intrinsic Wi^ ya distribution at 
1.9 < z < 2.7 and 2.7 < z < 3.4 is simply an equally- 
weighted combination of the distributions measured in 
these two redshift ranges. 11 In the bottom middle panel, 
the expected and observed distributions are not signifi- 
cantly different, so it is plausible that the intrinsic distri- 
bution of WLya for redshift 1.9 < z < 2.7 galaxies resem- 
bles an equally-weighted combination of the distributions 
measured for BX galaxies and LBGs. However, as the 
top middle panel indicates, such an intrinsic distribution 
overpredicts the number of galaxies at 2.7 < z < 3.4 with 
Lya in absorption. Decreasing the fraction of galaxies 

11 Since the comoving number densities of galaxies at redshifts 
1.9 < z < 2.7 and 2.7 < z < 3.4 are similar, it is reasonable to ap- 
proximate the combined WL ya distribution as an equally weighted 
sum of the measured distributions in these two redshift ranges. 



with absorption in the intrinsic distribution (e.g., by as- 
suming some non-equally weighted combination of Wi, ya 
distributions at low and high redshifts) may result in a 
better match for the observed z ~ 3 distribution, but 
would lose agreement with the observed z ~ 2 distribu- 
tion. 

Finally, if we assume that the intrinsic W\, ya distribu- 
tions at 1.9 < z < 2.7 and 2.7 < z < 3.4 are identical to 
those we actually measure in these two redshift ranges, 
then the expected distributions are very close to what is 
actually measured (right panels of Figure [8]). In fact, the 
expected distributions are only marginally different from 
the assumed intrinsic distributions, even if the assump- 
tions are erroneous (compare red and green rectangles in 
left and middle panels of Figure [H a Kolmogorv-Smirnov 
(KS) test indicates a > 50% probability that the intrinsic 
and expected distributions are drawn from the same pop- 
ulations). These observations suggest that the BX and 
LBG color selection criteria do not alter significantly the 
parent W^ ya distribution of galaxies. We already dis- 
cussed in § 13.31 why this must be the case for BX galax- 
ies, since there is a redshift range covered by BX selection 
where the U n GlZ colors are unaffected by Lya; we have 
just shown it to be true for LBGs also. 

In all cases shown in Figure [H we have quantified the 
disparity in the observed and expected W^ ya distribu- 
tions by computing the statistic £ as follows. We gener- 
ated 10000 realizations of the expected Wi, ya distribution 
(in the same way as we did to compute the uncertainties 
in the expected distribution; see above). We then per- 
formed a KS test to determine the probability (pks) that 
each of these realizations are drawn from the same distri- 
bution as the observed Wi, ya distribution. The quantity 
£ is then defined as the ratio of the number of realiza- 
tions where pks < 0.5 to the total number of realizations 
(10000). Low values of £ indicate that the expected and 
observed distributions are less likely to have been drawn 
from the same parent distribution. The values of £ are 
given in each panel of Figure [8l and support our con- 
clusion that the color criteria do not significantly alter 
the intrinsic Wi^ ya distributions. In summary, we find 
evidence that the fraction of emission-line galaxies (/20) 
appears to increase with redshift (Table |4|) and that such 
a trend is most likely not due to selection bias, as demon- 
strated by the differences in the expected and observed 
W^Lya distributions for galaxies at lower and higher red- 
shift (Figure El). 

4.2. E(B - V) Distributions 

A useful by-product of the maximum-likelihood 
method (§ I3.4.2|3.4~3]) is the distribution of galaxy spec- 
tral shapes, parameterized by the color excess E{B — V), 
corrected for incompleteness. Figure [S] shows the best- 
fit E(B — V) distributions, compared with the observed 
distributions, for galaxies at redshifts 1.9 < z < 2.7 and 
2.7 < z < 3.4 (al so tabulated in Table O. Using the V c g 
method (§ I3.4.1[) results in E(B — V) distributions that 
are within 10% of the observed distributions and there- 
fore deviate significantly from our best-fit distributions 
at z ~ 2 and z ~ 3. The analysis indicates that the true 
E(B — V) distributions are slightly bluer, on average, 
than observed. Table [5] lists the mean and dispersion 
of E(B — V) for the observed and maximum-likelihood 
distributions for galaxies at redshifts 1.9 < z < 2.7 and 
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Fig. 9. — Comparison of best-fit and observed E(B — V) distributions for galaxies at redshifts 1.9 < z < 2.7 (left panel) and 2.7 < z < 3.4 
(right panel). Histogram s denote the observed distributions, computed from the G — 1Z colors, and assuming a constant star formation 
model attenuated by the ICalzetti et a l, (2000) law. The width of each bar reflects the Poisson error in the corresponding bin. The red 
lines and yellow shaded regions indicate the mean and 1 a errors on the maximum-likelihood (best-fit) distributions. The blue dashed 
lines indicate the distribution without correcting for Lya perturbation to the observed colors. Dashed and solid vertical lines, respectively, 
denote the average E(B — V) for the observed and best-fit distributions. The E(B — V) distribution data are summarized in Table [5] 

2.7 < z < 3.4. 

The measured E(B—V) distribution for the BX sample 
is expected to be slightly biased toward redder spectral 
shapes than the intrinsic values because our photomet- 
ric method makes the colors appear slightly redder than 
they really are — and thus E{ B — V ) is redder — par- 
ticularly for fainter galaxies (§ I3.4.2p . In addition, the 
presence of Lya absorption in a galaxy's spectrum will, 
depending on the redshift, cause the G-band magnitude 
to appear fainter than the true broadband magnitude 
(corrected for line effects), such that E(B — V) will be 
overestimated. This latter effect can be visualized for 
1.9 < z < 2.7 galaxies in the left panel of Figure [9] the 
distribution uncorrected for the effects of Lya (dashed 
blue line) is systematically redder than the corrected dis- 
tribution (solid red line), owing to the presence of Lyct 
absorption among, and the low /20 value of, the major- 
ity of galaxies at 1.9 < z < 2.7 (e.g., Figure El El) . The 
systematic effects induced by Lya are less apparent in 
the E(B — V) distribution for 2.7 < z < 3.4 galaxies, 
primarily because the LBG selection window spans a re- 
gion of color space that is significantly larger than the 
typical color change induced by Lya perturbations. 

Before proceeding with a discussion of E(B — V) as an 
indicator of dust reddening and the variation of E(B — V) 
with redshift and apparent magnitude, we remind the 
reader that we can only correct for the incompleteness 
of objects whose colors are such that they are scattered 
into the selection windows. In other words, there are 
undoubtedly galaxies at these redshifts that will never 
scatter into the BX/LBG selection windows, for exam- 
ple, those galaxies that are optically-faint either because 
they have little star formation or are very dusty star- 
bursts (e.g., DRGs and SMGs). Therefore, the E(B-V) 
for such dusty galaxies will not be reflected in the dis- 
tributions shown in Figure [5] Typically, such very dusty 
galaxies would have E(B — V) > 0.45, although a sig- 
nificant fraction also show bluer E(B — V) comparable 



to those of BX/LBGs (jChapman et all [2005). Because 
these dusty star-forming and quiescent galaxies are in 
large part optically-faint, not accounting for them in our 
analysis should minimally affect our E(B — V) distribu- 
tion for optically-bright galaxies (Figure [5]). Further, as 
we show in § 18.21 comparison of our UV LF with those 
derived from magnitude limits surveys suggests that we 
must be reasonably complete for UV-bright (1Z < 25.5) 
galaxies at z ~ 2 — 3. 

4.2.1. E(B — V) as a Proxy for Dust Reddeniny 

Up until now, we have been using E[B — V) (the rest- 
frame UV slope) to parameterize the range of spectral 
shapes observed among high redshift galaxies. A number 
of studies have shown that E(B — V) also has a physical 
interpretation: it correlates very well with the redden- 
ing, or dust obscuration, of most high redshift galax- 
ies (e.g.. ICalzetti et ,al.ll200Ct lAdelberger fc Steide]||2000t 
iReddv et all l2006bf ) . Here we define reddening as the 



attenuation of luminosity by dust which can be parame- 
terized, for example, by the quantity Lboi/^uv- Combin- 
ing Spitzer MIPS data for a sample of spectroscopically- 
confirmed redshift 1.5 < z < 2.6 galaxi es where K- 
corrections could be computed accurately, IReddv et al.l 
(2006b) showed that E(B — V) not only correlates with 
Lboi for galaxies with Lboi ^$ 10 12 3 L Q , but that the cor- 
relation is identical to that established f or local galaxies 
(ICalzetti et al.ll2000t iMeurer et al.lll999l ). 

The E(B — V) for relatively dust-free (or very young) 
galaxies is dominated by intrinsic variations in the SEDs 
of high redshift galaxies, and so E(B — V) is not a direct 
indicator of reddening for these galaxies (which is why we 
measure a non-negligible number density of galaxies with 
E{B — V) < when assuming a single SED). Further, 
there is a significant presence of very dusty galaxies with 
-^bol ?t 10 12 L Q that are optically-bright (TZ < 25.5) and 
satisfy the rest-UV color criteria, but have E{B — V) that 
severely underpredict their attenuations and bolometric 
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TABLE 5 
Normalized E(B — V) Distributions 



E{B - V) 


BX (Measured) 


1.9 < z < 2.7 (Max-Lik) 


LBG (Measured) 


2.7 < z < 3.4 (Max-Lik) 


-0.2 0.1 a 


< 0.01 


0.02 ± 0.02 


0.02 ±0.01 


0.03 ±0.01 


-0.1 — 0.0 a 


0.04 ± 0.01 


0.10 ± 0.02 


0.08 ±0.01 


0.15 ±0.03 


0.0 — 0.1 


0.26 ±0.01 


0.32 ± 0.03 


0.26 ± 0.02 


0.30 ±0.05 


0.1 — 0.2 


0.47 ±0.02 


0.34 ± 0.02 


0.35 ± 0.02 


0.35 ±0.06 


0.2 — 0.3 


0.21 ± 0.01 


0.19 ±0.03 


0.21 ±0.01 


0.13 ±0.04 


0.3 — 0.4 


0.03 ±0.01 


0.03 ±0.02 


0.06 ±0.01 


0.03 ±0.02 


0.4 — 0.5 


< 0.01 


0.01 ± 0.01 


0.01 ±0.01 


0.01 ±0.01 


(E(B - V)) 


0.15 ±0.07 


0.12 ± 0.12 


0.14 ±0.07 


0.10 ±0.09 



a We measure a non-negligible number of galaxies with E(B — V) < since the E(B — V) of dust-free and/or 
very young galaxies is dominated by intrinsic variations in the SED (see text). 



luminosities (Re ddv et all l2006bf ). Finally, we remind 
the reader that we cannot account for the E(B — V) 
of objects that have a zero probability of being scat- 
tered into our sample. Nonetheless, our completeness- 
corrected estimates of the E(B — V) distributions suggest 
an average attenuation between dust-corrected and un- 
corrected UV luminosity, L^y/Luv, of ~ 4 — 5. This 
is similar to the value measured f rom (1) stacked X- 
ray data for BX galaxi es and LBGs (|Nandra et al.ll2002t 
iReddv fc Steidell 12004 . and (2) MIPS lum inosities and 
dust-corrected UV and Ha luminosities (|Reddv et al.1 
I2006bb lErb et al.ll2006bf). It is also the same value advo- 
cated bv lSteidel et all (|1999D in correcting observed UV 
luminosities for dust extinction among z ~ 3 LBGs. 

4.2.2. Comparison of Reddening Distributions with Redshift 

Remarkably, we find very little evolution in the red- 
dening distribution between redshifts 1.9 < z < 3.4, de- 
spite the roughly 730 Myr timespan between the mean 
redshifts for the low ((z) = 2.30) and high ((z) = 3.05) 
redshift samples, as shown in Figure [TO] It is not surpris- 
ing that the two distributions should span a similar range 
of E(B — V) since the BX criteria were designed to se- 
lected galaxies with a similar range of spectral properties 
as LBGs. Even so, the incompleteness corrections mod- 
ulate two very different observed E(B — V) distributions 
for the lower and higher redshift samples to the point 
where they are virtually identical. The difference in the 
fraction of large W\, ya emission systems between the two 
samples (§ 14. ip , and the possibility that such large Wh ya 
emission systems could be young and relatively dust-free 
galaxies (§ EJ) , does little to modulate the overall redden- 
ing distributions since such galaxies constitute a small 
fraction of the overall population. 

The similarity in E(B — V) and dust attenuation be- 
tween z ~ 2 and z ~ 3 galaxies agrees with the stacked 
X-ray studies of BXs and LBGs (jReddv fc Steidell 120041: 
iNandra et aLll2002f ). As we discuss in § the lack of 
evolution in E(B — V) implies that the extinction prop- 
erties of bright (1Z < 25.5) star-forming galaxies are not 
changing significantly between redshift z ~ 3 and z ~ 2, 
unlike the situ ation at lower (z < 2) and higher [z > 3 ) 
redshifts (e.g.. lReddv et alJl2006M iBouwens et al.ll2006fl . 

4.2.3. Reddening Distribution as a Function of Rest-Frame 

UV Magnitude 

Before turning to a discussion of the LF, we must 
first determine whether the best-fit E(B — V) distribu- 
tion shows any systematic changes as a function of rest- 
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Fig. 10. — Comparison of maximum-likelihood E(B — V) distri- 
butions for galaxies at redshifts 1.9 < z < 2.7 and 2.7 < z < 3.4. 



frame UV magnitude, since such changes can, in prin- 
ciple, affect the shape and normalization of the LF. As 
a first test, we restricted the maximum-likelihood analy- 
sis (§ 13.4. 213. 43)) to particular magnitudes in the range 
22.0 < 1Z < 25.5, and we did not find any significant 
trend in the E(B — V) distribution as a function of mag- 
nitude to 1Z = 25.5. 

Because E(B — V) can be used as a proxy for the 
reddening of galaxies (§ 14.2. ip . we have further inves- 
tigated whether the reddening distribution varies as a 
function of observed apparent magnitude by exploiting 
the multi- wavelength data in several survey fields. To 
accomplish this, we relied on our interpretation of the 
Spitzer MIPS data for a sample of BX-selected galaxies 
in the GOODS-N field; these data give us an independent 
probe of the dust emission in z ~ 2 galaxies. Figure QT] 
shows the dust obscuration factors, parameterized as 
Lboi/L\jy, where Lboi = ^m + iuv (infrared plus UV lu- 
minosity), 12 as a function of observed optical magnitude, 
from the MIPS ana lysis of the GOODS-North field by 
IReddv et all (|2006bD . The open red circles indicate rest- 
UV-selected objects at 1.5 < z < 2.6, most of which are 
BX galaxies, detected at 24 /j,m, and the large pentagon 
and crosses denote the average stack and distribution 
in 1Z magnitude, respectively, for galaxies undetected 

12 Here we define Lm as the total luminosity between 8 and 
1000 jim. 
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Fig. 11. — Distribution of attenuation factors, parameterized as 
Lb i/Luvi as inferred from Spitzer MIPS data, as a function of 
apparent optical magnitude 1Z for rest-UV-selected galaxies with 
redshifts 1.5 < z < 2.6. Also indicated is the stacked average 
for 48 galaxies undetected at 24 /im (large blue pentagon) and 
unconfused with brighter sources; the distribution in TZ magnitude 
for a larger sample of 73 galaxies undetected at 24 fim is shown by 
the arbitrarily normalized crosses. Total stacked 24 (im and X-ray 
averages are indicated by the solid red circles and green points, 
respectively, that include all galaxies. 



at 24 (tm. The total 24 /im stacked averages includ- 
ing both detected and undetected galaxies at 24 are 
shown by the solid red circles. Similarly, the Lboi/^uv 
inferred from X-ray stacked avera ges (computed in man- 
ner similar to that presented in iReddv fe Steidel 120041 
where Lboi is determined from the X-ray flux) including 
all galaxies, irr espective of direct dete ction in the Chan- 
dra 2 Ms data (| Alexander et al.ll2003t ) are shown by the 
green points. While there is some evidence that the dis- 
persion in attenuation factor increases towards fainter 
magnitudes, as evidenced by the larger spread of 24 //m 
detection galaxies and as would be expected if optically- 
faint galaxies have contribution from both heavily dust- 
obscured objects as well as those with intrinsically low 
star formation rates, the results of Figure [11] suggest that 
the average extinction correction (based on the stacked 
points) is approximately constant over the range in 1Z 
magnitude consi dered here. 13 These r e sults confirm the 
trends noted by lAdelberger fe Steidell (|2000l ). who used 
local templates to deduce that the observed UV lumi- 
nosities of galaxies at redshifts z = 0, z ~ 1, and z ~ 3, 
are inse nsitive to dust obscuration, Lboi/£uv (e.g., Fig- 
ure 17 of lAdelberger fe S tcidcl 20001). We have confirmed 
this trend explicitly at redshifts 1.5 < z < 2.6. The ob- 
served (unobscured) UV luminosity (i.e., the emergent 
luminosity after attenuation by dust) to 1Z = 25.5 will 
also be insensitive to bolometric luminosity since dust 
obscur ation is tightly corr elated with bolometric lumi- 
nosity (|Reddv et al.ll2006bT) . While the attenuation fac- 
tors and bolometric luminosities of z ~ 2 — 3 galaxies are 
insensitive to the unobscured UV luminosity, at least to 
TZ = 25.5, there is a very strong dependence of the dust- 



13 We note that Rcddy et al. (2006b) excluded objects from their 
analysis that were d irectly detected in the Chandra 2 Ms data in 
the GOODS-N field HAlexander et al.ll20"03fi of which almost all are 
AGN. 



corrected UV luminosity (or IR or bolome tric luminos- 
ity) on the attenuation factors of galaxies (jReddv et al.1 
I2006bt lAdelberger fc Steidel l2000t see also § GO}. For 
the purposes of the present analysis, we will assume 
that the reddening distribution of galaxies is constant 
to 1Z = 25.5. We return to the issue of how a varying 
reddening distribution affects our calculation of the total 
luminosity density in § 16.41 

5. RESULTS: UV LUMINOSITY FUNCTIONS 

5.1. Preferred LFs 

To provide the closest match between rest-frame wave- 
lengths, and thus avoid cosmological JT-corrections, we 
used 7?.-band as a tracer of rest-frame 1700 A emission 
at the mean redshift of the LBG sample, (z) ~ 3.05. 
Similarly, we used a "composite" magnitude between G 
and 7\L-band (man) as a tracer of rest- frame 1700 A at 
the mean redshift of the BX sample, (z) ~ 2.30, where 
rriQTi is simply the magnitude corresponding to the av- 
erage of the G and 1Z fluxes. Absolute magnitudes were 
computed using the standard relation: 

Mab(1700A) = m - 51og(d L /10pc) + 2.51og(l + 4(4) 

where Mab(1700A) is the absolute magnitude at rest- 
frame 1700 A, d^ is the luminosity distance, and m is the 
apparent magnitude at 7?.-band at z ~ 3 or at the com- 
posite G7<!.-band at z ~ 2. We have made the reasonable 
assumption that the SED JT-correction is approximately 
zero for the average rest-UV SED of BX-selected galax- 
ies after a star formation age of 100 Myr for the typical 
reddening (E(B — V) ~ 0.15) of galaxies in our sample. 

The maximum- likelihood rest-frame 1700 A luminos- 
ity functions for z ~ 2 and z ~ 3 galaxies are shown in 
Figure [121 and listed in TableGH The LFs were computed 
by using the entire photometric sample and holding the 
best-fit E(B — V) distribution (as determined from the 
spectroscopic sample; Figure G3) fixed. The extension of 
the spectroscopically-determined E(B — V) distribution 
to the photometric sample is a reasonable approximation 
given that (a) the spectroscopic and photometric samples 
are likely to have the same redshift distribution (§ I2.6P 
and (b) the E(B — V) distribution remains unchanged 
as a function of 1Z magnitude to 1Z = 25.5 (§ 14.2. 3jl . By 
nature of the maximum-likelihood method, our LF com- 
putation includes corrections for the systematic effects of 
photometric bias and Lya perturbations. Errors in the 
luminosity functions reflect both Poisson counting statis- 
tics and field-to-field variations; the latter are accounted 
for by examining the dispersion in the LF as a function of 
magnitude f or each of the field s of the survey (see § 15. 4|) . 
The best-fit ISchechten ()1976f ) function and parameters 
for the z ~ 2 LF are also indicated. Uncertainties in the 
faint-end slope a, characteristic absolute magnitude M* 
(or characteristic luminosity L*), and characteristic num- 
ber density cf>* are estimated by simulating many realiza- 
tions of the LF as allowed by the errors (assuming the 
errors follow a normal distribution), fitting a Schechtcr 
function to each of these realizations, and then determin- 
ing the dispersion in measured values for a, Al*, and <f>* 
for these realizations. 

Based on integrating our maximum-likelihood LFs, the 
fraction of star-forming galaxies with redshifts 1.9 < 
z < 2.7 and M AB (1700A) < -19.33 (i.e., 1Z = 25.5 
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TABLE 6 

Rest-Frame UV Luminosity Functions of 1.9 < 



< 3.4 Galaxies 



Redshift Range M AB (1700A) (xlCT 



h 3 
rl 0.7 



Mpc" 



ag 1 ) 



1.9 < z < 2.7 



Fig. 12. — Rest-frame UV luminosity functions at z ~ 2 (solid circles) and z ~ 3 (solid and op en squares fo r groun d-based observations 
and HST, respectively) computed in our analysis, compared with z ~ 4 results (triangles) from Stcidcl ct al. (1999). All data have been 
recast to the same cosmology used throughout this paper. Also indicated are the best-fit Schcchterl 1 119761) functions for the z ~ 2 (blue 
line) and 2 ~ 3 (red line) LFs. No shift in normalization was applied to the LFs. Confidence contours demonstrate the degeneracy between 
a and M* for the z ~ 2 fit, as shown in the inset. The red cross denotes a and M* for z ~ 3 galaxies. 

with redshifts 2.7 < z < 3.4 and Mab(1700A) < -20.02 
(7£ = 25.5 at z = 3.05) that have colors that satisfy the 
LBG criteria is « 47%. Note that some galaxies escap- 
ing LBG selection will be scattered into the BX window, 
and vice versa. Also, some galaxies that are intrinsically 
fainter (or brighter) than 1Z = 25.5 will be scattered into 
(or out of) the BX and LBG samples due to photometric 
error. The total fraction of galaxies with 1.9 < z < 3.4 
and 7Z < 25.5 that satisfy either the BX or LBG criteria 
is 0.55. 

In the following sections, we examine various as- 
pects of the luminosity functions derived here, including 
the differences in the LF derived using the V e g versus 
maximum-likelihood method, the significance (or lack 
thereof) of the Schechter parameters, and field-to-field 
variations. We conclude this section by examining how 
photometric redshifts can introduce non-trivial biases in 
the computation of the LF. 



2.7 < z < 3.4 



-22.83 - 


-22.33 


0.012 ±0.008 


-22.33 — 


-21.83 


0.05 ±0.01 


-21.83 - 


-21.33 


0.17 ±0.04 


-21.33 - 


-20.83 


0.53 ±0.08 


-20.83 - 


-20.33 


1.10 ±0.07 


-20.33 - 


-19.83 


2.13 ±0.07 


-19.83 - 


-19.33 


3.66 ± 0.02 


-23.02 — 


-22.52 


0.0031 ±0.0027 


-22.52 — 


-22.02 


0.025 ± 0.007 


-22.02 - 


-21.52 


0.09 ±0.01 


-21.52 — 


-21.02 


0.27 ±0.03 


-21.02 — 


-20.52 


0.60 ±0.03 


-20.52 - 


-20.02 


1.16 ±0.03 



at z = 2.3) that have colors that satisfy BX criteria is 
w 58%. Similarly, the fraction of star-forming galaxies 
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5.2. Comparison of the VeS an d Maximum-Likelihood 

Methods 

The maximum-likelihood technique was used to de- 
rive the LFs presented here. However, because many 
published LFs are derived using the less accurate V c g 
method, particularly for dropout samples at high red- 
shift, it is useful to determine how close (or how far) 
such determinations are from reality by comparing with 
our maximum-likelihood value. 
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Fig. 13. — Comparison of the maximum-likelihood LF at z ~ 
2 with the V e ff determinati ons from our work (open circles) and 
Sawicki & Thompson (2006) (open triangles). 

Figure [T3] compares the LFs at z ~ 2 computed using 
the V e s f§ !3.4.ip and maximum-likelihood (§ I3.4.2[) meth- 
ods, along wit h a comparison of the T4g d etermination 
at z ~ 2 from lSawicki fe Thompson] (|2006l) . The figure 
clearly demonstrates the systematic bias at both bright 
and faint magnitudes of the V e s LF with respect to the 
maximum-likelihood value. These biases are particularly 
apparent for criteria that target a narrow range in color 
space, such as the BX criteria, where photometric scatter 
or perturbations due to Lya can be as large or larger than 
the width of the color selection windows (e.g., Figure QJ 
[5]). For example, we would have inferred a significantly 
shallower faint-end slope of a = —1.21 ± 0.15 had we 
relied on the LF derived from the V c s method. 14 For 
the LBG criteria, we find little difference in the V e s and 
maximum-likelihood determinations of the LF, mostly 
due to the fact that the LBG selection window covers a 
larger area of color space and outlier objects (with col- 
ors placing them near the selection boundaries) for which 
the bias is the largest will make up a significantly smaller 
fraction of the sample. 

5.3. Schechter Parameters 

The spectroscopic sample allows us to accurately con- 
strain the LF, taking into account sample completeness, 

14 Note that our V e ff determ ination is slightly different from that 
of Sawicki & Thompson ( 2006 ) , despite the use of the exact same 
filter set and color criteria between the two studies, since our V c f{ 
determination includes the effects of Lyo line perturbations to the 
rest-UV colors and incorporates the maximum-likelihood E(B — V) 
distribution in the LF calculation. 



interloper fraction, and line perturbations, for galaxies 
with 7Z < 25.5. It is brighter than this limit that we 
consider our LF to be most robust. The results for z ^ 2 
galaxies fainter than 1Z = 25.5 is less certain given that 
our determination of the z ~ 2 faint-end slope relies on a 
spectroscopic sample that extends only ~ 4 times fainter 
than the characteristic luminosity of z ~ 2 galaxies. For- 
mally, we find a faint-end slope at 1.9 < z < 2.7 of 
a = —1.84 ± 0.11, with a characteristic magnitude of 
MX-a dJOOA) = -20 . 97 ± .23. 

The ISteidel et ail (| 19991 ) analysis of the z ~ 3 LF 
included [/-dropout galaxies in HDF-N where the red- 
sh ift distribution w as modeled using the color criteria 
of lDickinsonl (| 19981 ) and assuming the range of intrinsic 
spectr al shapes of LBGs found by lAdelberger fc Steidell 
(2000) (open squares in Figure ITS]) . Based on the com- 
bined I£eck_S£eeteoscopic and HDF-N [/-dropout sam- 
ples, ISteidel et al.1 (1999) found a steep faint-end slope 
a = —1.60 ±0.13. Further refinement of t he inc omplete- 
ness corrections by Adclbcr ger fc Steidell (|2000l ) resulted 
in a combined fit to the ground-based spectroscopic and 
HDF [/-dropout samples of a = -1.57 ± 0.11. Fitting 
only the ground-based (spectroscopically determined) 
points from our analysis at z ~ 3 yields a = — 1.85±0.22 
and MX B (1700A) = -21.21 ± 0.16. Combining our 
data with the HDF-N [/-drop points, and excluding the 
faintest HDF point that may suffer from incompleteness 
(|Steidel et alj|1999l ). results in a fit with a = -1.57±0.11 
and M£ B (1700A) = -20.84 ± 0.12. Not surprisingly, 
the se values are in excellent agr eement with those found 
bv lAdelberger fc Steidell (|2000t) . primarily because the 
same faint (HDF) data are used to determine the faint- 
end slope. The best-fit Schechter function at z ~ 3 is 
also indicated in Figure [T2J 

While a and M* are useful in parameterizing the gen- 
eral shape of the LF, we caution against over-interpreting 
their validity when accounting for faint galaxies that are 
beyond current spectroscopic capabilities. The absence 
of spectroscopic constraints on the asymptotic faint-end 
slope and a less than exponential fall-off of bright sources 
both conspire to make a steeper (i.e., more negative). 
However, these parameters are useful in describing a lo- 
cal approximation to data points that are not far from 
L*. For convenience, the best-fit parameter values from 
our analysis are listed in Table [7] In the subsequent anal- 
ysis, we will assume a = —1.6 for the faint-end slope of 
the UV LF at z ~ 2 — 3. Note that if the steeper faint- 
end slopes inferred from our shallower ground-based data 
(a ~ —1.85) accurately reflect reality, then this will serve 
to increase the total UV luminosity density of galaxies 
with SFRs between 0.1 and 1000 M Q yr" 1 by - 20% and 
~ 50% at z ~ 2 and z ~ 3, respectively, relative to the 
values obtained with a = —1.6. 

5.4. Field-to- Field Variations 

Access to multiple uncorrelated fields allows us to 
judge the effects of large scale structure on the derived 
LF. The dispersion in normalization between the lumi- 
nosity function in bins of 1Z derived in individual fields 
is a strong function of 1Z, as illustrated in Figure [T3] for 
the z ~ 2 sample. The points in this figure are deter- 
mined using the following steps. First, we computed 
the maximum-likelihood LF in each of the 14 fields of 
the z ~ 2 survey. The dispersion in LF values from 
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TABLE 7 

Best-fit Schechter Parameters for UV LFs of 1.9 < z < 3.4 Galaxies 





Redshift Range 


a 


MX B (1700A) 


<t>* 




x 2 


1.9 < 


z < 2.7 (ground-based) 


-1.84 ±0.11 


-20.97 ±0.23 


(1.74 ±0.63) x 10~ 


3 


0.73 


1.9 < 


z < 2.7 (ground-based) 


-1.60 (fixed) 


-20.60 ±0.08 


(3.31 ±0.22) x 10~ 


-3 


5.81 


2.7 < 


z < 3.4 (ground-based) 


-1.85 ±0.22 


-21.21 ±0.16 


(1.02 ±0.52) x 10~ 


-3 


0.18 


2.7 < 


z < 3.4 (ground+space) 


-1.57 ±0.11 


-20.84 ±0.12 


(1.66 ±0.63) x 10- 


-3 


5.72 



field-to-field, within a given magnitude bin, are deter- 
mined by weighting the LF values by the field size such 
that LF determinations from larger fields are given more 
weight than LF determinations from smaller fields. The 
fractional dispersion in normalization is then defined as 
the ratio of the dispersion of these weighted values and 
the weighted mean value of the LF in each bin. This 
fractional dispersion in normalization is much larger at 
the bright-end for 1Z < 22.5 and decreases significantly 
for galaxies with fainter magnitudes. This trend results 
from statistical fluctuations at the bright end due to 
the smaller number of galaxies and the fact that the 
clustering c orrelation function is a strong function of 
magnitude (Adel berger et all [2005a) . We further note 
that at least 4 of the 14 z ~ 2 survey fields show sig- 
nificant reds hift-sp ace over-densities (e.g., HS1700 field; 
ISteidel et al.l 12005). The effect of such over-densities on 
the derived LF will of course depend on the redshift of 
the over-densities with respect to the BX selection func- 
tion. An over-density at z = 2.8 is unlikely to affect the 
LF derived for 1.9 < z < 2.7 galaxies to the same extent 
as an over-density at z = 2.3 (placing it in the middle 
of the BX selection function) . One option when working 
in single fields is to use the available spectroscopy and 
known selection function to model the effects of such 
over-densities on the derived LFs, or use Monte Carlo 
simulations to estimate uncertainties in the normaliza- 
tion of the derived LF (|Bouwens et al.l [20061. Because 
our analysis includes many uncorrelated fields (14 and 
29 for the z ~ 2 and z ~ 3 samples, respectively) spread 
throughout the sky, we assume that the average LFs are 
representative of 1.9 < z < 3.4 galaxies. Any remain- 
ing uncertainty in normalization of the average LF (i.e., 
the uncertainty reflected in the field-to-field fractional 
dispersion shown in Figure I14[) is added in quadrature 
with Poisson counting error (shown as open circles in 
Figure I14[) to determine the total error bars shown in 
Figure [T2j We remind the reader that the systematic 
effects of photometric bias and Lya perturbations are 
already reflected in the derived LFs. 

5.5. Effect of Photometric Redshifts 

In light of recent literature regarding photometric es- 
timates of the UV LFs at redshifts z ~ 2 — 3 (e.g., 
iGabasch et alll200l . it is worthwhile to briefly exam- 
ine how our derived LFs would change in the absence of 
spectroscopic information, instead relying on photomet- 
ric redshifts (z p hot)- We derived the photometric redshift 
error, defined as 

Az = Zphot z spoc ; (15) 

for a sample of 925 star-forming galaxies with spectro- 
scopic redshifts 1.4 < z spcc < 3.5 that lie in fields with 
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Fig. 14. — Fractional dispersion in normalization of the z ~ 2 
UV LF as a function of apparent magnitude, both in terms of field- 
to-ficld (solid circles) and Poisson (open circles) variations. 

sufficient multi-wavelength data to warrant SED analy- 
sis. This is by far the largest spectroscopic sample at 
these redshifts, and it enables us to investigate how Az 
varies as a function of both redshift and magnitude since, 
in principle, the error will depend both on the relative 
placement of spectral breaks across the photometric fil- 
ters (i.e., the redshift of the galaxy) and on the quality 
of the photometry and significance of the detection (i.e., 
the apparent magnitude of the galaxy). 

Photometri c redshifts were estima ted using the Hy- 
perZ code of iBolzonella et ail ((2000). We only consid- 
ered galaxies with detections in at least the following 
bands: G,1Z, and K s . At least half of the resulting 925 
objects also have detections in either the J-band and/or 
Spitzer IRAC bands. All but 52 of the 925 objects are 
detected at U n ; the remaining 52 are all "C" candidates 
(§ 12 . 3[) at redshifts z > 2.7. We considered a variety of 
star formation histories, reddening, and redshifts when 
fitting the data using HyperZ. Figure [15] compares the 
photometric and spectroscopic redshifts for the sample 
of 925 galaxies. The biases and dispersions in photo- 
metric redshifts for galaxies fainter and brighter than 
M* = —20.97 in different redshift ranges are listed in 
Table both including and excluding catastrophic out- 
liers with Zphot < 0.6. Even excluding outliers with 
Zphot < 0-6 results in significant redshift error dispersions 
of cr(Az) ~ 0.33 — 0.42. Further, in all cases over the red- 
shift ranges where we compute the LFs, 1.9 < z < 2.7 
and 2.7 < z < 3.4, we find significant photometric biases 
ranging from Az ~ 0.2 — 0.5, in the sense that z p hot is 
systematically under-estimated (i.e., luminosity is over- 
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Fig. 15. — (Left:) Comparison between the photometric redshifts derived using HyperZ (Bolzonella ct al. 2000) and spectroscopic redshifts 
from our ground-based survey for a sample of 925 star-forming galaxies. The solid line denotes the unity relationship, and dashed lines 
demarcate the region over which the UV LF is computed. (Right:) Comparison of the UV LF derived using the photometric redshift 
distribution in the left panel (open circles) with our spectroscopic determination (solid circles). 



TABLE 8 

Biases and Dispersions in Photometric Redshift Errors 
(Az = z phot — 2 spC c) for Star-Forming Galaxies at 1.4 < z < 3.5 a 



Rcdshift-Range 


All Galaxies 
M > -20.97 M < -20.97 


Excluding z phot 
M > -20.97 


< 0.6 Galaxies 
AI < -20.97 


1.4 < z < 1.9 
1.9 < z < 2.7 
2.7 < z < 3.5 


0.01 ±0.40 
-0.26 ±0.58 
-0.50 ±0.93 


0.07 ±0.48 
-0.28 ±0.62 
-0.43 ±0.86 


0.04 ±0.35 
-0.18 ±0.42 
-0.20 ±0.42 


0.07 ±0.33 
-0.14 ±0.40 
-0.21 ± 0.39 



Az 



estimated). 

A simulation was constructed to examine the effect of 
these biases and dis persions on the LF, sim ilar to the 
method presented in iMarchesini et all (|2007| ). but modi- 
fied to (a) allow for many realizations of the intrinsic LF 
and (b) account for photometric redshift errors using the 
empirical data of Figure [15] and Table [H To accomplish 
this, we first generated many realizations of the UV LF as 
allowed by the (normally-distributed) errors of the spec- 
troscopically determined LF. We then drew magnitudes 
randomly from a Schechter distribution determined by 
fitting a Schechter function to each of these realizations 
of the intrinsic LF. Since the probability of a galaxy lying 
at redshift z, 



1 



dV 

p[z) oc — — oc , 

dz (i + z y ^/n A + n M (i + z) 3 



(16) 



is roughly constant over the redshift interval 1.4 < z < 
3.5, we drew redshifts from a uniform distribution. The 
result is a list of simulated redshift and magnitude pairs, 
(-Zspeo M), for galaxies. A photometric redshift was as- 
signed to each galaxy by randomly drawing a redshift 
from the distribution of z p hot (in left panel of Figure I15|) 
within a box of width 5z spcc — 0.4 centered at z spcc , 
using the z p hot distribution for galaxies either brighter 
or fainter than M* — —20.97 depending on the magni- 



tude M of the simulated galaxy. The absolute magni- 
tude of the galaxy was then recomputed assuming the 
photometric redshift. We then reconstructed the LF at 
1.9 < z < 2.7 for each realization assuming the photo- 
metric redshifts. The average LF from these many re- 
alizations is indicated by the open circles in the right 
panel of Figure [151 Because there are more galaxies 
scattered into, rather than out of, the redshift range 
1.9 < z < 2.7, the net result is that we would have 
over-estimated the intrinsic LF had we relied on photo- 
metric redshifts. While the difference in the photometric 
and spectroscopic LFs on the faint-end is small, it be- 
comes quite significant for galaxies brighter than M* . 
This systematic difference arises from the fact that the 
change in absolute magnitude (AM) for a fixed Az and 
apparent magnitude will be larger for galaxies scattered 
from low to high redshift than for galaxies scattered from 
high to low redshift. For example, a galaxy at redshift 
z S poc = 1.6 scattering to redshift z p hot = 2.0, implying 
\Az\ — 0.4, results in AM « 0.44. However, a galaxy at 
redshift ,z sp cc = 3.0 scattering to redshift z p hot = 2.6 (i.e., 
the same \Az\ as above) results in AM w 0.26. The net 
effect is that the bright-end of the LF is systematically 
inflated with respect to the faint-end. 

There are three further issues to note. First, it has 
become common in the literature to estimate photomet- 



22 



Reddy et al. 



ric redshift errors independent of fitting the stellar pop- 
ulations of galaxies by simply shifting prescribed fixed 
templates until a best-fit redshift is reached. Redshift 
errors derived in this manner will underestimate the true 
error in redshift obtained by marginalizing over the un- 
certainties of fitting those templates to the broadband 
photometry. Second, the simulation performed here ben- 
efited from the a priori knowledge that all the galaxies 
truly lie at the correct redshifts 1.9 < z < 2.7. Photo- 
metric redshift scatter (e.g., Figure [T5|) will generally be 
larger than that presented here since there will undoubt- 
edly be some very low redshift galaxies (z < 1.4) that 
are scattered into the range 1.9 < z < 2.7. Third, we 
remind the reader that the photometric redshift errors 
derived here are for optically-bright (7Z < 25.5) objects 
with spectroscopic redshifts. It is likely that the photo- 
metric redshift errors will be larger than assumed here 
for very faint galaxies where the photometric uncertain- 
ties may be larger. This, in turn, may bias the faint-end 
slope more severely than reflected in our simulations. We 
stress that the results of the photometric redshift simula- 
tion presented here (Figure [TS"]) are unique to our sample. 
As a result, the biases in the LF may be different for sur- 
veys that incorporate a different number of photometric 
filters with differing photometric data quality, although 
the photometric r edshift accuracy foun d her e is similar to 
that pr esented in lShaplev et al.l (|2005f ) and lReddv et al.l 
(2006b) using more (different) bands. In any case, this 
section illustrates how photometric redshifts can induce 
non-trivial biases in the LF. 

5.6. UV LF Summary 

To summarize, this section has focused on our measure- 
ments of the UV LF at z ~ 2 and z ~ 3. Our method 
for computing the LFs takes into account a number of 
systematic effects including contamination from low red- 
shift interlopers and AGN, Lya line perturbations to the 
observed colors of galaxies, and photometric scatter. A 
large number of independent fields allows us to control 
for sample variance. Further, spectroscopic redshifts en- 
able us to precisely correct for the effect of IGM opacity 
on the rest-UV colors. Our method for computing the 
LFs uses a maximum-likelihood technique to account for 
the systematic scattering of galaxies in parameter (e.g., 
luminosity, reddening, and redshift) space. Given this 
detailed treatment, we consider the UV LF at z ~ 2 and 
z ~ 3 derived here to be the most robust measurements 
yet. Comparison of our UV LF with the (corrected) de- 
termination from a magnitude limited sample suggests 
that our determination must be reasonably complete for 
galaxies with 1Z < 25.5 (see §12]). 

In the following sections, we will discuss how we can 
combine our determinations of the UV LFs at z ~ 2 — 3 
with what we know about the extinction properties of 
high redshift galaxies to infer LFs at other wavelengths. 
We will primarily focus on inferences of the IR and bolo- 
metric LFs at z ~ 2 — 3, but also present Ha LFs at 
similar redshifts, the latter of which may be useful for 
current and future emission line studies. 

6. RESULTS: REST-FRAME 8 (M, INFRARED, AND 
BOLOMETRIC LUMINOSITY FUNCTIONS 

As suggested in the previous analysis, correcting the 
rest-UV LF for the effects of dust extinction is a key 



component in recovering the star formation rate density. 
Aside from our knowledge of the E(B — V) distribution 
at high redshift (§ 14. 2[) . extensive multi- wavelength data 
enable us to independently determine the extinction cor- 
rections relevant for the same sample of spectroscopically 
confirmed galaxies that we used to compute the UV LF. 

Before the advent of panchromatic galaxy surveys, 
it was common to simply apply an avera ge correction 
for ex tinction, typically a factor of 4 — 5 (Steid eTet al.l 
1999). Subsequently, extensive multi- wavelength data 
have placed our extinction corrections on a much more 
solid footing. For i nstance, initia l X-ray and radio stack- 
ing a nalyses (e.g., iNandra et al.l l200"2l : iReddy fe Steidell 
l2004t ) indicated that high redshift UV-selected popula- 
tions with 1Z < 25.5 have average obscuration factors 
(Lir/Luv) aroun d 4 — 5, supporting t he average correc- 
tion advocated bv lSteidel et alj (|1999h . Further progress 
was made by taking advantage of the unique sensitivity 
of the Spitzer MIPS instrument, allowing us to directly 
detect for the fir st time the dust emi ssion fr om L* galax- 
ies at z > 1.5 (|Reddv et al.l l2006bf ). The IReddy et all 
(2006bj) analysis confirmed the average trends established 
by previous X-ray stacking studies, and further demon- 
strated that moderate lumi nosity galaxies ( 10 11 < Lboi ^ 
10 12 - 3 ) at z ~ 2 follow the IMeurer et alj (|1999) attenu- 
ation law found for local UV-selected starburst galaxies. 
The importance of this analysis for the present study is 
that we can directly relate the E(B — V) distribution 
of most z ~ 2 galaxies (§ 14. 2|) with their distribution in 
obscuration, Lir/Luv- 

In this section, we present our constraints on the 8 /im, 
infrared, and bolometric luminosity functions of z ~ 2 — 3 
galaxies, as derived from our UV LFs and the known 
extinction properties of galaxies at these redshifts. We 
present IR LFs based on two different recipes (our E(B — 
V) distribution and the distribution of 24 /im fluxes) for 
evaluating the dust attenuation of galaxies. In § 16. 1[ we 
combine our UV LF with the E(B — V) distribution to 
infer the IR LF. In subsequent sections, we combine our 
UV LF with the observed 24 /im properties of galaxies 
to infer the IR LF. The two methods are compared in 
detail in § El 

6.1. Extinction- Corrected Measures of the Luminosity 

Function 

As a first step, we can use the IMeurer et al.1 {1999) 
relation to recover the dust-corrected LF. The method 
proceeded with the following steps: 

1. We first generated many realizations of the 
maximum-likelihood UV LF and E(B — V) distribution 
at z ~ 2, assuming normal LF and E(B — V) distribu- 
tion errors. We randomly chose an LF and E(B — V) 
distribution from these many realizations, to create an 
IF/E(B-V) pair, {£,£}. 

2. Because the LF and E(B — V) distributions do not 
change significantly over the redshift range 1.9 < z < 2.7 
and because the E(B — V) distribution is insensitive 
to absolute magnitude down to our spectroscopic limit 
(§ 14.2. 3|) . we can assume that the intrinsic redshift z, 
magnitude M, and reddening E(B — V) of a galaxy are 
independent variables. Redshifts were drawn randomly 
from a uniform distribution. Magnitudes were drawn 
from the range —23 < Af(1700A) < —15.5 according 
to a Schechter distribution that described the luminos- 
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ity function C from the {£, £ } pair. The faint limit of 
M (1700 A) — —15.5 corres ponds to an unob scured SFR 
of - 0.1 M© yr _1 using the lKennicuttl (|1998ft calibration. 
We drew galaxies down to this low limit of unobscured 
luminosity because such galaxies can be scattered to bins 
of higher luminosity after correcting for extinction. Sim- 
ilarly, E(B — V) values were drawn randomly from the 
E(B — V) distribution £ from the {£, £ } pair, excluding 
negative E(B — V) values that reflect unphysical redden- 
ing values. The result is a list of galaxies associated with 
a triplet (z,M ,E(B - V)). 

3. The rest- frame 1700 A specific luminosity of each 
galaxy is calculated as 

47Td T 2 



Li 



700 



(1 



.^Q-0.4(48.60+m 1700 ) 



(17) 



where g?l is the luminosity distance at redshift z and 
mnoo is the apparent magnitude of the galaxy with ab- 
solute magnitude M at redshift z (Eq. [14")) . We then 
calculate vL v at 1700 A to yield the UV luminosity. 
The E( B — V) for the gala xy is used in conjunction 
with the lCalzetti et all (2000) relation to derive the dust- 
corrected UV luminosity. 

4. To determine the IR luminosity corresponding to 
this dust-corrected UV luminosity, we assumed that the 
UV and IR emission are tied directly to the SFR of the 
galaxy. The IR luminosity is assumed to be the lumi- 
nosity which, when added to the unobscured UV lumi- 
nosity, yields the same SFR that would have been ob- 
tain ed from th e dust- corrected UV luminosity, assuming 
the IKennicuttl l|1998f ) relations. These IR luminosities 
are then perturbed by a normal distribution with sigma 
of 0.3 dex to account for the dispersion between dust- 
corrected UV and IR luminosity (or, alternatively, the 
dispersion between E(B — V) and IR luminosity; e.g., 
IMeurer et al"1ll999t ). 

5. These IR luminosities are then binned to produce 
an IR LF. This is the IR LF corresponding to the {£, £} 
pair selected in step (1). 

Steps (l)-(5) are repeated many times, each time ran- 
domly drawing different {£,£} pairs. Aside from uncer- 
tainties in the rest-frame UV faint-end slope, there are 
two other systematics that can bias our determination 
of the IR LF: (1) a change in the faint-end slope of the 
rest-frame UV LF and (2) a change in the attenuation of 
UV-faint galaxies. We now discuss these two systematic 
effects in detail. 

First, we must determine how changing the number 
density of such faint objects, determined by the faint-end 
slope a, affects the IR LF. In principle, we could simply 
fix a — —1.6 when fitting the different realizations of the 
UV LF at z ~ 2 and compare with the results obtained 
by allowing a to vary freely in the Schechter fits to the 
realizations. However, this method will cause us to un- 
derestimate the errors on the faint-end of the IR LF. To 
obtain a truer estimate at the faint-end, we allowed a to 
vary freely around a normal distribution with a mean of 
{a) = — 1.6 and standard deviation of <x(a) = 0.11, sim- 
ilar to the dispersion in a that we measure when fitting 
the UV LF (Figure [H and Table [7]). For simplicity, we 
assume that a varies according to (a) = —1.6 ± 0.11 at 
z ~ 2 in the subsequent discussion. 

A second systematic effect that can bias the deter- 
mination of the IR LF is the distribution of extinction 



among TZ > 25.5 galaxies. For the calculation of the 
IR LF from E(B — V) we considered two cases. In the 
first case, we assume that the E(B — V) distribution 
is constant to arbitrarily faint rest-UV magnitudes. In 
the second case, we assume that the E(B — V) distri- 
bution is constant to TZ = 25.5 (e.g., Figure [IT] and 
§ I4.2.3|) . but then suddenly changes to have a mean 
(E(B — V)) = 0.04 (with same dispersion) for galax- 
ies fainter than TZ = 25.5 (we refer to this second case 
as a discontinuous E(B — V) distribution). We have as- 
sumed this value of E(B — V) = 0.04 because it is similar 
to the E(B - V) observed for very faint (< 0.1 L*) UV- 
selected galaxies inferred from dropout samples at higher 
rcdshifts (Bouwens et al. 2007, submitted). Because 
TZ = 25.5 is an arbitrary limit dictated by efficient spec- 
troscopic followup, it is highly unlikely that the E(B — V) 
distribution will suddenly change fainter than this limit. 
Rather, the distribution is likely to gradually fall towards 
lower E(B — V), or bluer rest-frame continuum spectral 
slopes (/3), proceeding to fainter galaxies, assuming that 
such fainter galaxies have lower star formation rates and 
are less dusty than 1Z < 25.5 galaxies. Therefore, the 
true E(B — V) distribution will very likely fall between 
the two extremes assumed above. We will return to this 
point shortly. 15 




Fig. 16. — Infrared luminosity fun c tions at z ~ 2 and z ~ 3, cal - 
culated assuming the IMeurer et all (fT999l ') and IKennicuttl p998h 
relations to convert unobscured UV luminosity and E(B — V) to 
infrared luminosities. The width of the shaded regions reflect the 
uncertainty in the rest-frame UV-slope and the attenuation distri- 
bution for 1Z > 25.5 galaxies (see text). 

For now, our IR LFs estimated from the E(B — V) 
distributions are shown in Figure [TH] and tabulated in 
Table ® The uncertainty in the LFs include uncertainty 
in the rest-frame UV faint-end slope and the uncertainty 
in the E(B — V) distribution for 1Z > 25.5 galaxies. The 
upper limit of each LF corresponds to the first case where 
E(B — V) is held constant. The lower limit of each LF 

15 We make the reasonable assumption that very dusty 
ULIRGs with faint UV luminosities make up a small fraction 
of the total number of UV-faint galaxies. Assuming other- 
wise would imply a significantly larger number of IR luminous 
galaxies than are presently observed in shallo w IR surveys (e.g., 
IPerez- Gonzalcz ct al. 2005: ICaputi et alj|2007h . 
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TABLE 9 

8 flM, IR, AND BOLOMETRIC LUMINOSITY FUNCTIONS OF 1.9 < 2 < 3.4 GALAXIES 3, 



1.9 < 2 < 2.7 2.7 < 2 < 3.4 (Predicted) 

(h 3 , 7 Mpc -3 decade -1 ) (h 3 , 7 Mpc -3 decade -1 ) 
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12.25 


- 12.50 


(1.10 ±0.34) 


X 


10- 4 


(6.55 ±2.37) 


x 10~ 5 


log L bo i b 


9.50 


- 9.75 


(3.69 ±1.78) 


X 


10 -2 


(2.57 ± 1.60) 


x 10~ 2 




9.75 - 


- 10.00 


(3.23 ± 1.27) 


X 


10~ 2 


(2.06 ± 1.15) 


x 10~ 2 




10.00 


- 10.25 


(2.47 ±0.96) 


X 


10~ 2 


(1.49 ±0.82) 


x 10~ 2 




10.25 


- 10.50 


(1.61 ±0.83) 


X 


10 -2 


(9.53 ±5.88) 


x 10~ 3 




10.50 


- 10.75 


(1.04 ±0.57) 


X 


10 -2 


(6.00 ±3.61) 


x 10~ 3 




10.75 


- 11.00 


(6.89 ±3.41) 


X 


10 -3 


(3.90 ±1.86) 


x 10" 3 




11.00 


- 11.25 


(4.25 ± 1.84) 


X 


10 -3 


(2.49 ± 1.05) 


x 10~ 3 




11.25 


- 11.50 


(2.51 ±0.80) 


X 


10 -3 


(1.48 ±0.45) 


x 10" 3 




11.50 


- 11.75 


(1.37 ±0.42) 


X 


10 -3 


(7.81 ± 1.98) 


x 10- 4 




11.75 


- 12.00 


(6.67 ± 1.32) 


X 


10 -4 


(4.05 ±0.74) 


x IO" 4 




12.00 


- 12.25 


(3.12 ±0.65) 


X 


10 -4 


(1.85 ±0.52) 


x IO" 4 




12.25 


- 12.50 


(1.17 ±0.38) 


X 


10- 4 


(7.67 ±2.12) 


x 10~ 5 



a Errors include systematic uncertainty in attenuation distribution for 1Z > 25.5 galaxies, 

as described in the text. 

b The values listed in this table are derived assuming the Caputi et al. (2007) calibration 
between vL v (S^ira) and Ljr. 



corresponds to the second case where E(B — V) suddenly 
decreases to have a mean of (E(B — V)) = 0.04 for galax- 
ies fainter than 1Z = 25.5. In general, the systematic un- 
certainties related to a changing attenuation distribution 
for UV-faint galaxies will dominate the uncertainties in 
the faint-end slope (§ 16. 3[) . 

6.2. Distribution of Dust Attenuation Factors 

As alluded to above, the distribution of rest-frame 
5 — 8.5 /im luminosities (Xs-s.s^m) of z ~ 2 galaxies 
observed by Spitzer/MIPS can be used to assess the in- 
frared luminosity function independent of any assump- 
tion regarding the relationship between rest-frame UV 
slope and extinction, as per the previous discussion. To 
this end, we must quantify the distribution of dust at- 
tenuation among z ~ 2 galaxies. For the subsequent dis- 
cussion, we will define the mid-IR (.Amir), far-IR (Air), 
and bolometric attenuation (Aboi) factors as the ratio be- 
tween uL v (8fj,m), Lir, and Lboi, respectively, an d £1700- 
Following the calibration of lReddv et al.l (|2006bf) . we can 



relate these A factors to each other: 
Air, « 12.9 Amir 

Aboi = Ll ^,;r" = A m + 1. (18) 

Note that these attenuation factors, A, are distinguished 
from the rest-frame UV attenuation factor which is the 
ratio between the dust-corrected and unobscured UV lu- 
minosities. 

The normal distribution of E(B — V) for galaxies at 
z ~ 2 — 3 (e.g., Figures [U [10]) implies that the atten- 
uation factors A will abide by a log-normal distribu- 
tion. We modeled the shape of the log A distribution 
by considering the measured log A of rest-UV-selected 
galaxi es with bolometric lu minosities Lboi < IO 12 3 L Q . 
From iReddv et all (|2006bh . (log Air) « 0.67, implying 
Lir/I/i7oo ~ 4.7, for the combined sample of 24 //m 
detected and undetected rest-UV-selected galaxies to 
1Z = 25.5. This mean attenuation implies the Gaussian 
fit shown in Figure 1171 compared with the distribution 
of log A for 24 /im detected galaxies. 



Multi- Wavelength LFs at 1.9 < z < 3.4 



25 




SMG 

(f(850yU,m)>5mJy): 



Fig. 17. — Distribution of measured log.4iR. for 24 fim de- 
tected galaxies with 7Z < 25.5 and luminosities Lbol < 10 12 ' 3 Lq, 
indicated by the red hashed histogram. The solid curve de- 
notes the Gaussian fit to the inferred distribution of all U n G7Z- 
selected galaxies with 1Z < 25.5, irrespective of 24 (im detection 
limit, and the vertical line indicates the mean of the distribution, 
(log Air) ~ 0.67. L og Air for bright SMGs from the analysis of 
Rcddv ct al. (2006b) is also shown. 

To construct a fair representation of the attenuation 
factors of high redshift galaxies, there is another is- 
sue that is pertinent. Namely, the distribution above 
does not take into account the non-negligible fraction of 
z ~ 2 — 3 galaxies that have attenuation factors much 
larger, on aver age, than those of typi c al galaxies at these 
redshifts (e.g.. IChapman et al.l l2003t iReddv et al.l l2005l: 
Ivan Dokkum et alj 120031 ). Virtually all of these galaxies 
have luminosities L bol > 10 123 L (jReddv et alj|2006bf) 
and « 50% of those that are also bright at submil- 
limeter wavelengths, /sso/xm ^ 5 mJy, also satisfy the 
BX/LBG criteria. Because our data are most sensi- 
tive to galaxies with luminosities Lboi J$ 10 12 ' 3 L© and 
because most galaxies with the largest attenuation fac- 
tors have I/boi ^ 10 123 L Q , we will only consider the 
^boi < 10 12 L Q regime when computing the IR LF. We 
combine our spectroscopically constrained estimate of 
the IR LF with higher luminosity data from the liter- 
ature in order to compute the total luminosity and star 
formation rate densities in § 18.31 18.41 and 18.51 



6.3. Attenuation of Rest- UV Faint Galaxies 

As discussed in § 16.11 the attenuation distribution of 
galaxies fainter than our 7?.-band limit for spectroscopy 
can affect significantly our inferences of the faint-end 
of the IR LFs. To remind the reader, in construct- 
ing the IR LFs that we would have inferred based on 
the UV continuum slope, we considered two cases. In 
the first, the distribution of E(B — V) is held fixed to 
the maximum-likelihood value (e.g., Figure [TO]) irrespec- 
tive of optical magnitude. In the second, we assumed 
the maximum-likelihood value for those galaxies with 
1Z < 25.5; for those fainter than this limit, we assumed 
a mean ( E(B — V)) = 0.04 ( corresponding to j3 ~ —2.0 
using the lCalzetti et aLll2000l relation). To construct the 
8 /im, IR, and bolometric LFs based on MIPS 24 /jm 
data, we assumed two cases similar to the ones consid- 
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25 26 27 
% (mag) 

Fig. 18. — Illustration of the two cases we consider for the at- 
tenuation distribution of galaxies as a function of rest-frame UV 
magnitude. In both cases, the height of the bars reflect a 1 a 
dispersion of 0.53 dex. 

ered above. In the first, we assume that all galaxies 
can be ascribed to the attenuation distribution shown 
in Figure [T71 In the second, we assume that the atten- 
uation distribution shifts to a mean of (log .Air) = 0, 
with the same dispersion as before, for galaxies with 
25.5 < 1Z < 27.5, then shifts again to a lower mean 
of (log Air) = —1 for galaxies fainter than 1Z — 27.5. 
The latter corresponds to an attenuation such that 90% 
of the bolometric luminosity of the galaxy emerges in 
the UV. Changing the distribution about the specified 
means does not significantly affect our conclusions. The 
two cases are illustrated in Figure [T8l 

It is very likely that the attenuation distribution does 
not remain constant to arbitrarily faint UV magnitude, 
as is assumed in case 1. Alternatively, because our 
1Z = 25.5 spectroscopic limit is arbitrary, we do not ex- 
pect the attenuation properties of galaxies fainter than 
this limit to drastically change. It is more likely that the 
attenuation distribution gradually shifts to lower mean 
values for galaxies fainter in the UV. Case 2 can therefore 
be considered a reasonable lower extreme to the gradient 
of attenuation as a function of rest-UV magnitude. The 
true IR LF derived from a realistic attenuation distribu- 
tion will likely lie between the IR LFs derived assuming 
case 1 and case 2. 

Is there any reason to believe that the attenuation 
has a much steeper gradient than "case 2" (i.e., bot- 
tom panel of Figure [TS]) ? Future deep stacking analy- 
ses and deep spectroscopic campaigns should resolve this 
question. For now we point out that much more ex- 
treme cases for the attenuation of UV- faint galaxies (i.e., 
such that they are even less obscured than what we have 
considered in case 2) would result in total IR LD and 
SFRD estimates that are comparable, if not lower, than 
estimates based on samples of the brightest (1Z < 25.5 
and K s (Vega) < 22) objects at these redshifts, most 
of which have luminosities comparable to LIRGs and 
ULIRGs (|Reddv et al.ll2005l ). Since the total IR LD and 
SFRD must be at least as large as that contributed by 
LIRGs and ULIRGs, the attenuation distribution of sub- 
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Fig. 19. — 8 /im and IR luminosity functions at redshift 1.9 < z < 2.7, compared with predictions in the higher redshift range 
2.7 < z < 3.4. The right panel also demonstrates the similarity in the IR LF derived from the rest-frame UV slope (Figure I16H with 
that derived using t he attenuation factors calculated from Spitzer/MlPS observations of 1.5 < 2 < 2.6 galaxies in the GOODS-N field 
flReddv et al.N2006bl) 

ULIRG galaxies cannot be much lower than what we have 

considered here. We will return to this point in § 18.51 TABLE 10 

Contributions of the IR LD at z ~ 2 



6.4. Rest-Frame 8 /im and IR Luminosity Functions 

Combining our UV LFs with the attenuation distri- 
bution derived from MIPS 24 /jm observations (using 
the same method described in § 16. ip results in estimates 
of the IR LFs. Given the amount of focus in using 
MIPS 24 /zm observati ons to probe star- f o rming popu- 
lations at z ~ 2 (e.g., ICarjuti et all l2007t IReddv et all 
l2006ttlPanovich et al J 1 200^ 7 it is useful to derive a rest- 
frame 8 /im luminosity function; we did this using the 
8 /zm attenuation f actors observed for z ~ 2 galaxies 
(jReddv et al.ll2006bl ). 16 We then use the relationship be- 
tween Am and .Amir (Eq. [TH]) to infer the IR luminosity 
function. 

Our inference of the 8 /im and IR LFs at 1.9 < z < 2.7 
at faint and moderate luminosities (Lir < 10 12 L©) are 
shown in Figure [T9l and listed in Tabled For later com- 
parison, we have assume d the relation b etwee n 8 /im and 
IR luminosity given by ICaputi et al.l (|2007h . The up- 
per and lower limits of the shaded regions in the figure 
correspond to the two different cases of attenuation dis- 
tributions discussed in the previous section. We cannot 
directly measure the rest-frame mid-IR luminosities of 
galaxies at 2.7 < z < 3.4, but we show the predicted 
8 /im and IR LFs at these redshifts assuming (a) the 
same attenuation distribution and (b) the same relation- 
ship between mid-IR and total IR luminosities found for 
z ~ 2 galaxies. Figure [19] demonstrates that the IR LF 
we would have inferred at 10 9 5 < Lir S 10 12 L© from 
the rest- frame UV slope, or E(B — V), is consistent with 
the one inferred from the MlPS-determined attenuation 
factors of these galaxies. This similarity reflects the sig- 
nificant correlation between E(P> — V) and attenuation 

16 For ease of comparison with previous literature, we express 
the mid-IR attenuation factor »4ir in terms of uL ^(8tim) rather 
than Z/5_ 8.5 Mm as was used in IReddv et al.1 l!2006bP ). The relation- 
ship between the two for the typical mid-IR SED of star-forming 
galaxies is uL^(8fj,m) Si 0.75L5_8.5 M m- 



Lir 



log IR LD 



10 9 - 10 10 L 

10 10 _ 10 n Lq 

10 11 — 10 12 L 
> 10 12 L 

Total (6 x 10 s < Lir < oo) 1 



8.03 ±0.17 (9%) 
8.45 ± 0.09 (23%) 
8.68 ± 0.06 (39%) 
8.48 ± 0.32 (25%) a 
9.09 ± 0.08 



a From ICaputi et al.1 S2007t) . 

b The lower limit of Lir = 6 X 10 s L roughly corre- 
sponds to an SFR of 0.1 M yr _1 . 

for galaxies with moderate luminosities, and such galax- 
ies are typical of the redshift z ~ 2 — 3 population (i.e., 
with luminosities corresponding to ~ L*). 

Because our data are most sensitive to galaxies with 
Lir ^ 10 12 L , we must incorporate direct measure- 
ments of the IR LF for high luminosity objects, such as 
those from Spitzer mid-IR surveys. There are several 
published values of the IR LF for ULIRGs ; here we as- 
sume the most recent determination from ICaputi et al.l 
(|2007f h Table [TU] lists the contribution of galaxies in 
different luminosity ranges to the IR LD, where we 
take the total LD to be that of galaxies with Lir > 
6 x 10 8 L . This limit corresponds to galaxies with SFRs 
of ~ 0.1 M Q yr . The total LD changes negligibly by 
integrating to zero luminosity. Table [TU] shows that - 
despite the large systematic uncertainties at the faint-end 
induced by variations in the attenuation distribution — 
a significant fraction of the IR LD at z ~ 2 arises from 
galaxies with sub- ULIRG luminosities. We will return 
this point in § 18.41 

6.5. Bolometric Luminosity Functions 

Finally, to gain an accurate picture of the distri- 
bution of total energetics of star- forming galaxies, we 
must consider the combined contribution from unob- 
scured (UV) luminosity and obscured (IR) luminosity. 
While the bolometric luminosity should closely follow 
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the infrared luminosit y for luminous galax ies (Xboi ~ 
Lm. > 3 x 10 11 L Q ), iReddv et all (|2006bl ) show that 
such an assumption is no longer valid for galaxies with 
Lho\ ^ 3 x 10 11 L (at z ~ 2), given the very tight 
correlation between dust obscuration and bolometric lu- 
minosity. For example, a 10 11 Lq galaxy at z ~ 2 will on 
average have half of its total luminosity emerging at UV 
wavelengths; similarly, a 10 10 ' 5 Lq galaxy at z ~ 2 will 
on average have 84% of its bolometric luminosity emerg- 
ing at UV wavelengths. 17 Figure [5TJ] shows the bolometric 
luminosity functions of star-forming galaxies at redshifts 
1.9 < z < 3.4, and values are listed in Table [9l The 
bolometric LF is larger in all luminosity bins considered 
than the IR LFs given that objects will shift from lower to 
higher luminosity bins after accounting for the emergent 
UV luminosity of high redshift galaxies. We note that 
in computing our prediction for the bolometric LF at 
2.7 < z < 3.4, we have assumed the same distribution of 
attenuation factors that was found for z ~ 2 galaxies, as 
was done in computing the IR LFs. The bolometric LFs 
are presented here because they give a true picture as to 
the total energetic output of galaxies, irrespective of dust 
extinction or the fraction of unobscured luminosity. In 
section § 18.51 we will discuss what the spectroscopically- 
constrained bolometric LFs imply for the contribution of 
moderate luminosity (10 11 < iboi S 10 12 Lq) galaxies 
to the global luminosity density. 
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Fig. 20. — The bolometric luminosity function at redshift 
1.9 < z < 2.7, calculated using the sum of the UV (unobscured) 
and IR (obscured) luminosities of galaxies. Our prediction of the 
bolometric LF at 2.7 < z < 3.4 is also shown. 



17 While the tight correlation between Z/bol ar >d attenuation has 
been observed both locally and at high redshift, the normaliza- 
tion of the relationship increases at higher redshift. This means, 
for example, that galaxies at z ~ 2 are on average 10 times more 
luminous for a given dust obscuration (or are 10 times less dust ob - 
scured for a given -Lbol) than z = galaxies (Rcddy et al. 2006b). 
Hence, the fraction of total luminosity emerging in the UV is larger 
at higher redshift than it is locally for galaxies of a given bolomet- 
ric luminosity. Note that this observation is still consistent with 
the finding that dustier systems dominate t he luminosity density 
at z ~ 1 — 2 relative to the present-day fe.g,. lLe Floc'h et al .1120051 ; 
ITakeuchi et~aLi r2005). as we discuss in § 18.41 



7. RESULTS: Ha LUMINOSITY FUNCTION 

We briefly discuss our derivation of the Ha LFs here, 
as they may be useful for current and future high red- 
shift emission line studies. The key ingredient that al- 
lows us to convert our UV LF into an estimate of the 
Ha LF is the correlation between dust -corrected UV and 
Ha estimates of star formation rates. lErb et aTl (|2006bl ) 
found a significant (6.8 a) correlation with 0.3 dex scat- 
ter between the extinction corrected UV estimates and 
Ha es timates of the SFRs, assuming the ICalzetti et al.l 
(2000) relation, for a sample of 114 rest-frame UV se- 
lected galaxies at z ~ 2. We can invert the relationship 
between extinction-corrected SFR and Ha line emission 
in order to infer the Ha LF. The spectroscopic Ha obser- 
vations used to establish the correlation betwe en Ha and 
UV-det ermined SFRs are described in detail in lErb et al.l 
(|2006aH cT). 

7.1. Method 

The method used to estimate the Ha LF of z ~ 2 galax- 
ies is analogous to that presented in § 16.11 We generated 
many realizations of the LF and E[B — V) distributions 
and randomly selected magnitudes and E(B — V). To 
determine the Ha luminosity corresponding to this dust- 
corrected UV luminosity, we assumed that the UV and 
Ha emission are tied directly to the S FR of the g alax; 
where the SFR is calibrated using the iKennicuttJ (|199. 
relations. It is then easy to show that 



L HQ [crgs s !] « 1.77 x 10 13 Li 700 [< 



ergs s" 1 Hz _1 ](19) 



The resulting Ha luminosities are then perturbed by 
0.3 dex to account for the dispersion in the rela- 
tion between the dust-corrected UV and Ha estimates 
(|Erb et al.ll2006bl) . For consistency with previous deter- 
minations of the Ha LF at lower redshifts, it is useful 
to derive an Ha luminosity function uncorrected for ex- 
tinction at z ~ 2. To accomplish this, we assume that 
the E(B — V) value, which is derived from the rest-frame 
UV colo rs, reflects the ne bular reddening of the galaxy- 
see al so lErb et al.l (|2006bf )). Applying the lCalzetti et all 
2000) relation to the intrinsic Ha luminosity, and as- 
suming the same value of E(B — V), yields an estimate 
of the observed Ha luminosity. 18 Note that because the 
attenuation A\/E(B — V) is a factor of w 3 smaller at 
the wavelength of Ha, 6563 A, than at 1700 A, it is not 
correct to simply convert the observed UV luminosity 
to a SFR and then back to an observed Ha luminosity: 
one must take into account the differential extinction, for 
a given E(B — V), between these two wavelengths. The 
dust-corrected (intrinsic) and uncorrected (observed) Ha 
luminosities are then binned to produce a dust-corrected 
and observed Ha LF, respectively. We considered the 
same two cases for the E(B — V) distribution of UV- 
faint galaxies as in § 16.11 Our predictions for the uncor- 
rected and dust-corrected Ha LFs at z ~ 2 are shown in 
Figure I2T1 and listed in Table HU 

7.2. Predicted Ha LFs at z - 3 



18 lErb et all H2006bl ) show that assuming the same E(B — V) in 
dust-correcting the Ho esti mates, as opposed to a smaller nebular 
reddening as advocated by Calzetti ct al. (2000), results in better 
agreement with dust-corrected UV and stacked X-ray estimates. 
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Fig. 21. — Dust-corrected and observed (uncorrected for extinction) Ha LFs inferred for star-forming galaxies with redshifts 1.9 < z < 2.7 
(left) and 2.7 < 2 < 3.4 (right). The uncertainty of each LF, represented by the width of the shaded regions, is dictated by the u ncer tainty 
in (1) the E(B — V) distribution for 1Z. > 25.5 galaxies and (2) the rest-frame UV faint-end slope. Values are tabulated in Table [TT] 



TABLE 11 

Ha Luminosity Functions of 1.9 < z < 3.4 Galaxies' 1 





1.9 <z < 2.7 


2.7 < z < 3.4 (Predicted) 




Observed <j> Dust-Corrected <j> 


Observed <j> Dust-Corrected <j> 


log L(Hot/ergs s 1 ) 


(Mpc -3 decade -1 ) 


(Mpc -3 decade -1 ) 



41.00 


— 41.25 


(4.34 ± 1.66) x 10" 


2 


(4.44 ± 1.75) x 10" 


2 


(1.72 ±0.55) 


x 10" 


2 


(2.27 ± 1.09) 


x 10" 


2 


41.25 


- 41.50 


(3.07 ±1.29) x 10" 


-2 


(3.52 ± 1.71) x 10" 


2 


(1.33 ± 0.48) 


x 10" 


2 


(1.63 ±0.81) 


x 10" 


2 


41.50 


- 41.75 


(1.98 ±0.88) x 10" 


-2 


(2.51 ± 1.41) x 10" 


2 


(9.18 ±3.83) 


x 10" 


3 


(1.11 ±0.59) 


x 10" 


2 


41.75 


— 42.00 


(1.20 ±0.59) x 10- 


-2 


(1.67 ± 1.07) x 10- 


2 


(6.07 ±2.70) 


x 10" 


3 


(6.87 ±3.85) 


x 10" 


3 


42.00 


— 42.25 


(6.95 ±3.42) x 10- 


-3 


(1.05 ±0.75) x 10" 


2 


(3.94 ± 1.26) 


x 10" 


3 


(4.25 ±2.31) 


x 10" 


3 


42.25 


- 42.50 


(4.07 ± 1.62) x 10" 


-3 


(6.48 ±4.28) x 10" 


3 


(2.44 ± 0.82) 


x 10" 


3 


(2.70 ± 1.22) 


x 10" 


3 


42.50 


- 42.75 


(2.31 ±0.63) x 10- 


-3 


(3.93 ± 2.16) x 10- 


3 


(1.55 ±0.33) 


x 10" 


3 


(1.63 ±0.49) 


x 10" 


3 


42.75 


- 43.00 


(1.23 ±0.20) x 10" 


-3 


(2.35 ±0.99) x 10- 


3 


(7.63 ± 2.21) 


x 10" 


4 


(8.60 ± 1.94) 


x 10" 


4 


43.00 


- 43.25 


(4.80 ±0.66) x 10" 


-4 


(1.27 ±0.40) x 10" 


3 


(3.30 ± 0.93) 


x 10" 


4 


(3.91 ± 0.90) 


x 10" 


4 


43.25 


- 43.50 


(1.63 ±0.58) x 10" 


-4 


(6.22 ± 1.57) x 10" 


4 


(1.38 ±0.91) 


x 10" 


4 


(1.38 ±0.43) 


x 10" 


4 


43.50 


- 43.75 


(4.74 ± 1.71) x 10" 


-5 


(2.61 ±0.88) x 10" 


4 


(3.44 ± 2.23) 


x 10" 


5 


(5.13 ±2.41) 


x 10" 


5 


43.75 


— 44.00 


(1.34 ±0.58) x 10" 


-5 


(9.53 ±3.03) x 10" 


5 


(2.11 ± 1.89) 


x 10" 


5 


(1.24 ±1.01) 


x 10" 


5 



a Errors include the systematic uncertainty in the E(B — V) distribution for 1Z > 25.5 galaxies, as described in the 
text. 



The correlation between UV and Ha SFRs has only- 
been tested directly at redshifts 2.0 < z < 2.6, where 
the Ha line falls within the if-band, making it accessi- 
ble to near-IR spectr ographs, such as Keck II/NIRSPEC 
(jMcLean et al.|[l998l ). It is useful, nonetheless, to predict 
the form of the Ha LF at z ~ 3 assuming that the cor- 
relation between UV and Ha SFRs holds at these higher 
redshifts. This is a reasonable assumption to make given 
that, to 1Z = 25.5, the z ~ 2 and z ~ 3 samples host 
galaxies with a virtually identical range of E(B — V) 
fFigurefTD]) , and galaxies in these respective samples have 
similar average dust attenuation factors, represented as 
the ratio of the dust-corrected UV and unobscured UV 
luminosities, of ~ 4 — 5 based on X-ray and radio stack - 
ing analyses (|Nandra et al.ll2002tlReddv fc Steide]||2004D . 
Our predictions for the uncorrected and dust-corrected 
Ha LFs at z ~ 3, computed using the steps above and 
using the combined ground-based and HDF samples to 
generate the UV LF realizations, are shown in Figure [2T1 
and listed in Table [TT1 We briefly present a comparison 



of our Ha LFs and luminosity densities with others from 
the literature in § 18.3.21 

8. DISCUSSION 

We have presented the most robust estimates of the 
rest-frame UV LFs and moderate luminosity regime of 
the IR LFs of star forming galaxies at redshifts 1.9 < 
z < 3.4 (§ [7l [6j) . We have demonstrated how pho- 
tometric redshifts over this redshift range can introduce 
non-trivial biases in the LF, underscoring the need for 
spectroscopy where it is feasible (§ 15. 5| . Further, our 
extensive spectroscopy allows us to examine other sys- 
tematic effects, including Lya line perturbations to the 
intrinsic rest-frame UV colors of galaxies (§ I3.3l4.ljl . In 
the next section, we will examine more closely the trend 
between Lya emission and redshift and its dependence on 
the physical properties of galaxies. In § 18.21 we discuss re- 
cent results that have indicated an excess of bright galax- 
ies akj^JSjyvjii that o bser ved in the initial L BG studies 
of ISteidel etaU (| 19991 ) and iDickinsonl (|1998l ). Compari- 
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son of our UV, Ha, and IR LFs with those of previous 
studies, and the evolution in the LF and luminosity den- 
sity, are discussed in § 18.31 and 18.41 Integral constraints 
on the star formation rate density are presented in § 18.51 

8.1. Wj^ya Distribution as a Function of Redshift 

The analysis of § 14.11 demonstrated that the intrinsic 
W-Lyc distribution of galaxies at 1.9 < z < 3.4 was not 
significantly modulated by the BX and LBG color crite- 
ria and, furthermore, that the fraction of galaxies with 
W^Lya > 20 A (/20) was larger at higher redshifts. This 
trend was recognized by comparing /20 between (a) BX 
galaxies and LBGs and (b) BX galaxies at z < 2.48 and 
z > 2.48, with results summarized in Table 2] In the 
latter case, the fact that we see a trend in /20 with 
redshift even for galaxies selected using a single set of 
color criteria (BX) further strengthens our conclusions 
that /20 increases with redshift, irrespective of selection 
biases (see Figure |3|). 

To interpret this trend in a physical context, we as- 
sembled the stellar population parameters for galaxies 
with me asured Wlvc where SEP modeling was avai l- 
able from lShaplev et all (|2005f ) and lReddv et all (|2006af ). 
The resulting sample includes 139 galaxies, 14 with 
W^ya > 20 A. We used KS-tests to determine whether 
the SED parameters (star formation histories, ages, stel- 
lar masses, and star formation rates) for galaxies with 
Whya < 20 A are drawn from the same parent pop- 
ulation as those with Wh ya > 20 A. Doing this, we 
found no significant differences in the star formation 
histories, ages, stellar masses, and star formation rates 
of galaxies between these two samples. This result is 
not surprising given (1) the small sample size analyzed 
here, (2) the significant systematic degenerac i es between 
SEP parameters (e.g.. IShaplev et all 120051 : lErb et all 
I2006d: iPapovich etafl I2001D . and (3) the large uncer- 
tainty in the measured Wi^ ya for individual galaxies. 
Galaxies with W\, Ya > 50 A have an average age of 
~ 300 ± 300 Myr whereas those below this limit have 
an average age of ~ 460 ± 600 Myr. The difference in 
average age is not significant given the large dispersion 
in ages measured for the two samples. 

Nonetheless, several previous studies at z > 4 
suggest that Lya emitting galaxies are young (< 
50 Myr), low- metallicity systems with small stel- 
lar masses (e.g., IStanwav et al.l 120071: iPentericci et all 
20071: iPow-Hveelund et all l2007t iPirzkal et all 120061: 
Finkelstein et all 120071 : iLehnert fc Bremeri I2003D . par- 
ticularly in relation to galaxies w ithout Lya in emis- 
sion at the same redshifts (z ~ 6; IStanwav et al.1 120071 
iPow-Hvgelund et all 120071 ) . Results at z ~ 2 also in- 
dicate that low stellar mass and low metallicity galax- 
ies have significantly stronger Lya emission than high 
stellar mass and high metallicity galaxies (Er b et al.l 
2006al). Among UV - contin uum s elected sample s, both 
Dow-Hygelund et all (|2007l) and IStanwav et all (|2007l ) 
find a fraction of z ~ 6 galaxies wit h Wlvq 5, 25 A sim - 
ilar to the fraction found at z ~ 3 (jShaplev et all 120031 ). 
but the former find an excess of high equivalent width 
galaxies (WL yct > 100 A) compared to the z ~ 3 sam- 
ple. Further, the connection between Lya profiles and 
the physical properties of galaxies is well-known to be 
quite complicated, with a sensitivity to ionizing flux, dust 



obscuration, and v e locity of outflowing material (e.g., 
Taoken et all 120071 iReddy et all l20061i lHansen fc Ohl 
20061 : IShaplev et al.l [20031: lAdelbereer et al.ll2003lh De- 

spite these complications, the advantage of our large 
spectroscopic analysis is that we can very accurately 
quantify the trends between Lya emission and redshift 
(§ 14. ip . The presently small samples at z ~ 6 pre- 
vent us from determining whether the trends observed 
at 1.9 < z < 3.4 extend to higher redshift, but the gen- 
eral expectation is that if younger and less dusty galaxies 
preferentially show Lya in emission, then the frequency 
of such Lya emitting galaxies should increase with in- 
creasing redshift as th e average galaxy ag e (and average 
dust-to-gas ratios; see lReddv et al.ll2006bll decreases. 

8.2. VVDS- Inferred Excess of Bright 
M (17001) < -22.0 Galaxies at Redshifts 
2.7 < z < 3.4 

Recently, iLe Fevre et all ([2005) reported results from 
the /-band magnitude limited VIMOS VLT Deep Sur- 
vey (VVDS) that indicated a significant excess of bright 
(M(1700A) < -22.0) galaxies with redshi fts 2.7 < z < 
3.4 co mpared with the earlier results of ISteidel et al.1 
(Il999fl and those i nferre d here from our likelihood analy- 
sis. iPaltani et al.1 (|2007h further quantified this excess by 
casting it into the form of a luminosity function, which 
we reproduce in Figure l22l 19 

We suggest two reasons for the discrepancy between 
our LF and that of the VVDS. First, the frequency of ob- 
jects with redshifts outside the redshift range 2.7 < z < 
3.4 (i.e, the contamination fraction; / c ) for M(1700A) < 
—22.0 is sign i ficant l y larger than the valu e inferred by 
IPaltani et all (|2007l ): lLe Fevre et al.l(|2005l ). based on the 
statistics of our much larger spectroscopic sample. To il- 
lustra te this, we must first consider how IPaltani et al.l 
l|2007h weight their galaxies in their computation of the 
LF to account for f c . The VVDS redshifts used to com- 
pute the LF fall within 4 categories. Flag-1 and Flag- 
2 objects are considered to have the least secure red- 
shifts, whereas flag- 3 and flag-4 obje c ts are more secure 
(|Paltani et al.ll2007| ). IPaltani et all (|2007h then deter- 
mine a contamination fraction of / c ~ 0.54 for the 254 
flag- 1/2 sources and assume a value of f c = for the 
12 flag-3/4 sources that are used in their LF compu- 
tation. Weighting the fractions according to the num- 
ber of sources then yields a net contamination rate of 
[254 x 0.54 + 12 x 0]/[254 + 12] w 0.52. 20 

The actual contamination rate among VVDS objects 
must be larger than / c = 0.52 for several reasons. First, 
the BX and LBG criteria account for the U n GTZ col- 
ors of - 70 - 80% of VVDS obj ects claimed to lie 
at 2.7 < z < 3.4 (see Figure 5 oflFaltaniet al.l 120071 
and Figure 3 of ILe Fevre et al1l2005l where most of the 
objects lie in the same region of color space encom- 
passed by either the BX or LBG criteria). Yet, 77% 
of spectroscopically-confirmed candidates in the BX and 

19 Although the VVDS UV LF of lPaltani eFaXI (t200W) is com- 
puted in a slightly different redshift range, 3 < z < 4, from that 
considered for our z ~ 3 LF (2.7 < 2 < 3.4), the comparison be- 
tween the two is valid since we find little evolution in the number 
density of M(1700A) < —22.0 galaxies between z ~ 3 and z ~ 4. 

20 This gives the same result as the "photometric rejection" 
method, where 137 of 266 objects are lik ely to lie outside the red- 
shift range 2.7 < z < 3.4, as discussed in Paltani et al. ( 2 0071 ). 
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LBG samples with M(1700A) < -22.0 are low redshift 
interlopers (Table SJ). 21 

Second, VVDS objects that do not satisfy the BX and 
LBG color criteria because of their redder G — TZ col- 
ors lie in the same region of c olor space as low red- 
shift star-forming: galaxies (e.g.. lAdelberger et al.ll200"^ : 
iReddv et all [2005). Without additional K s -band data 
to exclude these low-z interlopers (e.g., BzK selection 
of iDaddi et alj l2004bf ). the contamination rate among 
these red G — TZ objects is likely to be at least as large 
as the rate among objects that are targeted by criteria 
specifically designed to selected galaxies at 2.7 < z < 3.4 
(e.g., the LBG color criteria). To quantify this further, 
we turned to the magnitude limited Team Keck Red- 
shift Survey database in the GOOD S-N field (TKRS; 
IWirth et all l200i ICowie et alJl200l . There are 2471 
TKRS sources with TZ < 24 (roughly corresponding to 
the VVDS magnitude limit) with matching U n GTZ pho- 
tometry in our catalog. Of these matches, there are 755 
sources that satisfy the following conditions: (a) spectro- 
scopically observed in the TKRS, (b) do not satisfy either 
the BX or LBG criteria, and (c) have U n — G > 0.4 and 
G — TZ < 1.8. These limits define the region in color space 
of the « 20% of VVDS objects that do not satisfy the BX 
or LBG criteria. 22 Of these 755 bright sources that were 
observed in TKRS, 581, or « 77%, are spectroscopically- 
confirmed to lie at redshifts z < 1.4. This is a strict 
lower limit to the contamination rate since there will be 
(a) some 1.4 < z < 2.7 galaxies that are missed sim- 
ply because the lines used for redshift identification are 
shifted out of the TKRS wavelength coverage and (b) 
some galaxies with z < 1.4 that are unidentified because 
of poor weather, bad reduction, and other reasons. 

Finally, the instrumental setup used by the VVDS to 
obtain spectroscopy, resulting in observed wavelengths 
of 5500 and 9500 A is highly non-optimal for select- 
ing LBGs, particularly because of the lack of coverage 
around the Lya line. In summary, our value of f c = 0.77 
is s ignificantly la r ger th an the value of f c = 0.52 claimed 
by iPaltani et al.1 (|2007t ). yet is based on our sample of 
285 secure spectroscopic redshifts versus only 12 secure 
redshifts and 254 unsecure redshifts of the VVDS sur- 
vey. For the reasons discussed above, the real contam- 
ination rate among VVDS objects is likely to be larger 
than f c = 0.77. Even so, assuming this fraction would 
conservatively lower the brightest VVDS LF points by 
a factor of at least two. We already alluded to in § 15.51 
how catastrophic redshift errors can artificially boost the 
bright-end of the LF. 

The excess of bright galaxies in the VVDS LF is also 
likely to be due in part to the presence of non star- 
forming galaxies at redshifts 2.7 < z < 3.4, in other 
words QSOs and other bright AGN. We have explicitly 
excluded AGN from our LF determination, as described 
in § 12.51 but this was not done for the VVDS analysis. 
Based on our spectroscopy, the AGN/QSO contamina- 
tion rate is a strong function of magnitude and is larger 

21 For the BX sample, most of the contamination at bright mag- 
nitudes arises from foreground galaxies. For the LBG sample, most 
of the contamination arises from stars. 

22 Note that although nominally the BX criteria are designed to 
select galaxies with z < 2.7, we include them in the discussion here 
since a significant fraction of VVDS objects claimed to lie at z > 3 
fall in the same region of color space as BX candidates. 



than ~ 60% in the brightest luminosity bin of the VVDS 
analysis (Tabled]). In fact, of the 14 spectroscopically- 
confirmed LBGs in our sample with M (1700A) < -22.5, 
11 (or roughly 80%) are QSOs or bright AGN. Applying 
the AGN/QSO contamination fractions (as determined 
from our spectroscopy) , will further reduce the brightest 
VVDS LF points by a factor of 2.5. 

10" 3 F— I ' ' ' 1 ' ' ' ' 1 ' ' ' ^ 




M AB (1700A)-51ogh 07 

Fig. 22. — Comparison of our z ~ 3 UV LF determination (circles 
for data and solid line for best-fit Schechter function) and that of 
the VVDS (open squares and dashed line). Correcting the VVDS 
points for (a) objects outside the redshift range 2.7 < z < 3.4 and 
(b) AGN/QSOs at redshifts 2.7 < z < 3.4, as determined from our 
spectroscopic sample, results in much better agreement between 
the two LFs, as indicated by the filled squares. 

To conclude our comparison, the two primary causes 
for the excess of bright galaxies inferred by the VVDS is 
(a) their underestimate of the fraction of objects that lie 
outside the redshift range 2.7 < z < 3.4 and (b) the frac- 
tion of bright AGN and QSOs at redshifts 2.7 < z < 3.4. 
Our spectroscopy, which is the most extensive at the red- 
shifts in question and was obtained using the most opti- 
mal instrumental setup to identify galaxies at 2.7 < z < 
3.4, allow us to very accurately quantify the magnitude 
of both sources of contamination, all within a combined 
survey area that is roughly twice as large as the VVDS 
field. Figure[22] demonstrates that applying the contami- 
nation fractions determined from our spectroscopy to the 
VVDS points (after factoring out the VVDS contamina- 
tion correction) results in a better agreement between 
the VVDS and our LF. Taking all these results into con- 
sideration, we find no convincing evidence for an excess 
of bright galaxies at 2.7 < z < 3.4. 

Finally, the agreement between our corrected estimate 
of the VVDS LF (derived from a magnitude-limited sur- 
vey) with the one from our likelihood analysis (derived 
from a color-selected survey), strongly suggests that our 
LF must be reasonably complete for UV-bright galaxies. 

8.3. Evolution of the UV and Ha Luminosity Functions 
and Densities 

8.3.1. UV LFs 

Figure [12] summarizes our determinations of the rest- 
fra me UV LFs at redshi fts 1.9 < z < 3.4, compared with 
the ISteidel et all (|1999f ) LF at z ~ 4. As our method 
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Fig . 23. — (Left:) Comparison of a few UV LFs from the litera ture. The local UV LF (derived from GALEX data) is shown (Wvd eFet al.l 
2005) along with the z ~ 1 UV LF from[Arnouts ct al. (2005), as well as higher redshift determinations from the B, V, and / dropout 
samples of Bouwens et al. (2007; submitted). (Right:) Unobsc ured (uncorrected fo r extincti on) UV luminosity density , integrated to a 
fixed luminosity of L ^ m = 0.04L*_ 3 , from the following sources: IWvder et all J2005T) at z ~ 0, Schiminovich ct al. (2005) at z ~ 0.3 - 1.0, 
Arnouts ct al. (2005) at z ~ 2.0 from a photometric redshift analysis of the HDF North and South fields, and Bouwens et al. (2007; 
submitted) at z ~ 4 — 6. Our determinations at z ~ 2 and z ~ 3 are shown by the large red pentagons. 



TABLE 12 

Summary of Total UV, Ha, and IR Luminosity Densities at 1.9 < z < 3.4 



UV LD a Ha LD IR LD b 

Redshift Range (ergs s" 1 Hz -1 Mpc -3 ) (ergs s _1 Mpc -3 ) (L@ Mpc -3 ) 

1.9 < z < 2.7 (3.63 ± 0.40) x 10 26 (4.41 ± 1.10) x 10 40 (1.23 ± 0.21) x 10 9 

2.7<z<3.4 (2.29 ± 0.34) x 10 26 (2.71 ± 0.57) x 10 4() (6.61 ± 1.24) x 10 8 



a Uncorrected for extinction and integrated to 0.04L*_ 3 . 

b Values assume the Caputi et al. (2007) calibration between vL v (&iim) and Lir. The IR 
LD at z ~ 2 includes the contribution from ULIRGs (Caputi ct al. 2007). 



of constraining the reddening and luminosity distribu- 
tions takes into account a number of systematic effects 
(e.g., contamination fraction particularly at the bright 
end of the LF, photometric bias and errors, Lya line 
perturbations to the observed colors 23 ) that have not 
been considered in previous analyses (e.g. JGabasch et ail 
12004 iLe Fevre et al.l I2005D or were only partially con- 
sidered flSteidel et al. 119991: lAdelberger 'fc Steidel l2000t 
iSawicki fc Thompsonll200l . we regard our LFs as the 
most robust determinations at z ~ 2 and z ~ 3 to 
ft = 25.5. 

Our analysis indicates that the rest-frame UV LF 
shows little evolution for galaxies brighter than M* (at 
z ~ 2) between redshifts 1.9 < z < 4.5: the number den- 
sity of galaxies brighter than M* — —20.97 appears to 
be constant over the ~ 1.3 Gyr timespan between z ~ 4 
and z ~ 2.3. This lack of evolution in the bright-end of 
the UV LFs does not specifically address how a galaxy 
of a particular luminosity will evolve. For example, the 
lack of evolution at the bright end (M AB (1700A) < -21) 
of the LF does not imply that there is a population of 
UV-bright galaxies that is unevolving. Rather, if galaxies 
follow an exponentially-declining star formation history, 
then UV-bright galaxies at z ~ 3 will become fainter in 
the UV by z ~ 2, but will not necessarily be absent from 
the z ~ 2 sample. A precipitously declining star forma- 
tion history may imply that some UV-bright galaxies at 

23 Our determinations of the LFs are insensitive to small changes 
in the assumed Wx,y a distributions, such as those caused by trends 
in Wijya with apparent magnitude and color. 



z ~ 3 will be too faint to be included in UV-selected 
samples at z ~ 2. In any case, the lack of evolution 
at the bright-end of the UV LF implies that whatever 
UV-bright galaxies at z ~ 3 — 4 fall out of UV-selected 
samples by z ~ 2 must be made up in number by younger 
galaxies, those that are merging and just "turning on" , 
and/or those that are caught in an active phase of star 
formation at z ~ 2. The net effect is that the number 
density of galaxies with (Mab(1700A) < —21) remains 
essentially constant. 

For galaxies fainter than M*, we do find evidence for 
a small evolution between z ~ 3 and z ~ 2: the num- 
ber density of galaxies fainter than M* = —20.97 is sys- 
tematically larger at z ~ 2 than at z ~ 3 to a lower 

luminosity limit of 11 = 25.5 or M AB (1700l) 19, 

although the result is not of great significance given the 
generally overlapping error bars on the points between 
z ~ 2 and z ~ 3. What is clear is that the number den- 
sity of -21 < M A b(1700A) < -19 galaxies at z - 2 is at 
least as large as the the corresponding number density at 
z ~ 3. To put these results in context, Figure [23] sum- 
marizes our results at z ~ 2 and z ~ 3 with a few results 
at higher and lower redshifts. 24 

Figure [221 demonstrates the evolution of the UV LF. 
The number density of bright galaxies increases from 
z ~ 6 to z ~ 4 causing a brightening of M* by 

24 The figure is not meant to be comprehensive with respect to 
all determinations of the UV LF, particularl y at z > 4 where there 
are some differences between studies (e.g ., Sawicki & Thompson 
[20061; IQuchi et aLll200ll; llwata et ani2Wt iBeckwith et al.11200611 . 
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Fig. 24. — (Left:) Comparison of our inferred dust -corrected Ha lumi n osity function at z ~ 2 and predicted one at z ~ 3, with 
the direct Ho LF determinations at lower redshift from [T resse et al. (2002); Tresse & Maddox (1998). (Righ t:) E x tinction-corrected Ha 
lumin os ity density: o p en squares are from Hopkins ( 2004: ) and inclu d e dete r minations from IGalleeo et al. (1995); Perez-Gonzale z" et alJ 
j2003) : I5ullivan et all J2000h; ITresse k, Maddoxl 119981) ; ITresse eFall (I2002T) ; IGlazebrook et al.l (11999): IHopkins et al.l (12000) ; lYan eTall 
(1999); Mo orwood et al.l 120001) . The point at z = 2.75 is the H/3 determination from lPettini et al.l (119981) . The large red pentagons denote 
values from this work. 



~ 0.6 — 0.7 mag. In contrast, the number of faint galax- 
ies appears to undergo less evolution between z ~ 6 and 
z ~ 4. To quantify this further, we have calculated the 
rest- frame 1700 A unobscured UV luminosity density 
(LD). The luminosity density and error are estimated 
by simulating many realizations of the UV LF consistent 
within the normally-distributed LF errors and evaluating 
the mean luminosity-weighted integral of the LFs from 
each of these realizations. The mean value and dispersion 
of these integrated values give the mean luminosity den- 
sity and error, and values are listed in Table [T2"l We have 
assumed a faint-end slope of a = —1.6 ± 0.11 (i.e., as in 
§ 16. 1|) at z ~ 2. For consistency, all the LFs are integrated 
to a luminosity limit of L\ lm — 0.04L* =3 . This corre- 
sponds to a luminosity of Lr lm = 3.9 x 10 27 ergs s _1 Hz -1 
at 1700 A, is equivalent to « 0.05Z/* =2 > an d corresponds 
to an unobscured SFR of ~ 0.5 M© yr _1 assuming the 
iKennicuttJ (|1998h relation. The right panel of Figure [23] 
summarizes the integrated (unobscured) UV LD as a 
function of redshift including several published values, 
showing the significant evolution between z ~ 2 and 
z = 0. 

While the observed evolution in the UV LF is hardly 
surprising, we have placed this evolution during the most 
active epoch of star formation (z ~ 2 — 4) on a secure 
footing with our extensive spectroscopic analysis and de- 
tailed completeness corrections. Our analysis covers a 
larger number of uncorrelated fields than what has typ- 
ically been considered in previous studies, thus enabling 
us to mitigate the effects of sample variance. We note 
that part of the evolution of the unobscured UV LF may 
be a result of extinction, as we will discuss shortly. 

8.3.2. Ha LFs 

Figure[M]compares our inference of the Ha LF at z ~ 2 
and our prediction at z ~ 3, with the direct deter mina- 
tions at lower reds hifts by ITresse fc Maddoxl (|1998l ) and 
ITresse et al.l ((2002) . The evolution of the dust-corrected 
Ha LF qualitatively mimics the evolution observed in the 
UV LF (Figure [53]) for redshifts z < 2 - 3. We see a fac- 
tor of two decline in the number density of moderately 
luminous galaxies with 10 41 < L(Ha) < 10 42 5 ergs s _1 
from z ~ 2 to z ~ 0.7. The decline of moderately lumi- 



nous galaxies is at least a factor of 4 — 10 between z ~ 2 
and z ~ 0.2. 

The systematic excess of z ~ 2 galaxies with respect to 
z ~ 3 for luminosities L(Ha) < 10 43 ergs s _1 is primarily 
a result of the fact that galaxies on the faint-end of the 
UV LF (where we observed the same small systematic 
excess; see Figure IT2"]) are scattered to correspondingly 
more luminous bins of Ha luminosity after correcting for 
extinction. The significance of the systematic excess at 
z ~ 2 with respect to z ~ 3 is hard to judge since, unlike 
the case with the UV LFs, there are no direct determina- 
tions of the Ha LFs at z > 2. What is certain is that the 
number (and luminosity) density of moderately Ha lumi- 
nous galaxies at z ~ 2 is at least as large as the number at 
z ~ 3. In contrast, the significance of the increased fre- 
quency of moderately luminous galaxies at z ~ 2 with 
respect to z ~ 0.7 is supported by the fact that the 
Ha luminosity density shows a decline from z ~ 2 to 
z = 0. To illust rate this, we ha ve compiled estimates of 
the Ha LD from lHopkinsI (|2004l ) in the right panel of Fig- 
ure [2H including our new inferences at z ~ 2 a nd z ~ 3 
(the la tter are listed in Table [T2"]) . Values from iHopkinsI 
(2004) have been reddening-corrected assuming a lumi- 
nosity dependent obscuration correction. Our values at 
z ~ 2 and z ~ 3 are computed from the dust-corrected 
Ha LFs. Quantitatively, the Ha LD per comoving Mpc 
decreases by a factor ~ 25 between z ~ 2 and the local 
value. While H/3 observation s have been perfor med for 
a handful of objects at z ~ 3 (]Pettini et al.|[l99^ ). larger 
samples are needed to directly constrain the Ha LF at 
z ~ 3. Since our Ha results on the SFRD are degenerate 
with those estimated from the UV and IR LFs, we will 
not discuss the Ha LFs any further. 

8.4. Average Extinction and IR and Bolometric LFs 
8.4.1. Evolution of Dust Obscuration with Redshift 



As mentioned in § 18.3. 1[ the evolution observed in 
the unobscured rest-frame UV LF between z ~ 2 and 
z = may be partly a result of systematic differ- 
ences in extinction with redshift. It is already known 
that galaxies of a fixed bolometric luminosity have an 
average attenuation factor at z ~ 2 that is ~ 8 — 
10 times smaller than the attenuation factors of lo- 



Multi- Wavelength LFs at 1.9 < z < 3.4 



33 



10 



12 _ 



z = 2 

(Reddy et al. 2006) 



© 



+ 

i — i 



o 



46%/ 12% 



10 



n _ 



10 



10 



10* 




z=l 
(Burgarella 
et al. 2007) 



z = 

(Buat et al. 2006) 



i 







i 



i i i 



0.1 



10 



100 



Lir/^uv 



Fig. 25. — Trend b etween bolomet ric luminos i ty an d dust attenuation based on the analyses of IReddv et al.l ((2006b ) at z ~ 2, 
IBurgarella et af] (I2007T I at z ~ 1, and Buat ct al. (2006) at z ~ 0. The thickness of the lines show schematically the contribution of 
galaxies in different ranges of bolometric luminosity to the bolometric luminosity density, with the rel ative fractions listed. Fractional 
contributions at z ~ 2 were determined from our estimate of the bolometric LF at this epoch (Figure |20| |. The contributions at z ~ and 
Z ~ 1 were taken from Lc Floc' h~et~aT] (I2005T) . 

very strong evolution in the LF between z ~ 2 and z ~ 
0. We find a similar evolution in the IR LFs and IR 
LDs, shown in Figures [26l and [27l respectively, where we 
summarize the comparison between the z ~ 2 and z ~ 3 
IR LFs and luminosity densities from our analysis with 
other published measurements. The important point is 
that the average extinction correction needed to recover 
bolometric LDs from UV LDs will depend both on the 
offset in the Lboi versus attenuation trends between z ~ 2 
and z = 0, as well as the relative numbers of galaxies in 
different luminosity ranges that contribute to the LD. 
For example, LIRGs at z ~ 2 are on average 8 — 10 times 
less dusty than LIRGs at z = 0, but there are many 
more LIRGs at z ~ 2 than at z ~ 0. The contribution of 
galaxies in different luminosity ranges to the bolometric 
LD at z ~ 2, z ~ 1, and z = are shown schematically 
in Figure [25l 

If the IR LF was unevolving between z ~ 2 and z = 0, 
then the fact that galaxies of a fixed Lboi are 8—10 
times less dusty at z ~ 2 than at z ~ would imply that 



cal galaxies (IReddv et alj|2006bt IBurgarella et all 12007 
Bua t et all 120071 ). Figure 1231 ada pted from Redd v et al.l 



l!u; 

(2006 b]) and the GALE X results of lBuat et all (|2006t ) and 
IBurgarella et al.1 (|2007l ). illustrates the offset between the 
z ~ 2 , z ~ 1 , and z = tre nds between Lboi and attenua- 
tion. IReddv et al.l (|2006bD interpreted this trend as a re- 
sult of the increasing extinction per unit SFR (or increas- 
ing dust-to-gas ratio) as galaxies age. The dependence 
of attenuation on bolometric luminosity at l ow redshifts 
has been discussed by man y authors (e.g.. iBuat et al.l 



rAfonso et al.l 120031 : 
Sullivan et aTlfe OOl: Hopkins et al. 2l)01;i,B. The analy - 



20051 Pcrcz-Gonzalcz e£iD l2003t 



ses of IReddv et al.l (|2006bD ; lAdelberger fc Steidell (|2000D 
and Figure l25l demonstrate that this dependence contin- 
ues unabated from z = 0toz~2 - 3. 

If the LF was unevolving between z ~ 2 and z = 0, 
then the offset shown in Figure [25] implies that, when 
integrating the UV (or Ha) LF to a fixed luminosity, 
the extinction correction will be larger at lower redshifts. 
However, as Figures [23] and [M] demonstrate, there is a 
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Fig. 26. — Comparison of our inference of the 8 fim (left) and IR LFs (right) at z ~ 2 and our prediction a t z ~ 3, with direct 
measurements from the literature: z ~ 1 and z ~ 2 points from|Caputi ct al. (2007) and Pcrcz-Gonzalcz ct al. (2005 j), 8 am LF at z ~ 0.2 
from H uang et al. (2007, submitted), and the local IR LF fr om the IRAS BGS s ample ISanders et al.H2003l) . The IP crcz-Gonzalcz ct al. 
(2005) points at z ~ 2 are significantly larger than those of Caputi ct al. (2007), primarily because the former exclude only the most 
extreme AGN from their analysis and adopt a conversion between rest-fra me mid-IR and IR luminosity that has be en shown to overpredict 
the bolometric luminosities of > 2 X 10 12 Lq galaxies by a factor of > 2 (Papovich ct al. 2007; Daddi ct al. 2007b; see text). 



the average extinction correction needed to recover the 
bolometric LD from the UV LD would be roughly 8 — 10 
times larger at z ~ than at z ~ 2. However, because 
of the evolution in number density, the average correc- 
tion at z ~ is not 8 — 10 tim e s larg er than the value 
at z ~ 2. S chiminovich et ail (|2005[ ) find an average 
attenuation factor for UV-selected samples at z < 1 of 
(^boi/iuv) ~ 7, which is only 1.6 times larger than the 
value of 4.5 found for UV-selected galaxies at z ~ 2 — 3. It 
is interesting to note that even taking into account the 
evolution in number density, the average correction at 
z ~ is still larger than the correction at z ~ 2, despite 
a greater fraction of the LD in dustier galaxies at high 
rcdshift. This result may be partly due to the redshift 
evolution in L D oi versus attenuation trend (Figure I25p , 
but the larger correction at z ~ may also be needed to 
account for IR luminous galaxies that altogether escape 
UV-selection at z ~ 0. 

8.4.2. Evolution of the Dust-Obscured and Bolometric 
Luminosity Densities 

There are several important conclusions to draw from 
our analysis of the IR luminosity density. First, the un- 
obscured UV LD drops by a factor ~ 23 between z ~ 2 
and z ~ 0, whereas the IR LD drops by a factor of 
~ 14 between z ~ 2 and the present-day. This differ- 
ence can be partly accounted for by an evolution in the 
average extinction correction, as discussed in the previ- 
ous section. The independent determinations of the av- 
erage extinction cor rections at z ~ 2 ((ibol/^uv) ~ 4.5; 
Reddv et all l2006bt) and at z ~ ({L ho i/L V y) ~ 7; 
Schiminovich et all l2005f ) result in an evolution of the 
dust-corrected UV LD that is in remarkable agreement 
with the evolution of the IR LD between z ~ 2 and z = 
(see SlOl. 

Second, the IR LD at z ~ 2 appears to be at least 
as large as the value at z ~ 3, a result consistent with 
that of the UV LD analysis. This should not come as 
a surprise since the same incompleteness-corrected rest- 



frame UV LF was used to infer the IR LF. Nonetheless, 
our spectroscopic analysis puts the constraints on the Ha 
and IR LD on a more secure footing. In particular, our 
analysis provides the first spectroscopic constraints on 
the sub-ULIRG regime of the IR LFs at z > 2. 

Third, the constraints on the IR LD between 1.9 < 
z < 3.4 yield critical information on the relative con- 
tribution of galaxies to the IR LD as a function of lu- 
minosity. Figure [26] hints that previous studies of the 
IR LF at z ~ 2 underpredicted the number (and lu- 
minosity) density of sub-ULIRGs with Lir < 10 12 L Q 
(i. e., compare our IR LF poin ts wi th the extrapolation 
of iPerez-Gonzalez et al.l (|2005| ) and lCaputi et all (|2007t ) 
for galaxies with Lm < 10 12 L Q in Figure [26]) . When 
integrating the IR LF to account for all galaxies, we 
find a total IR LD of - 1.2 x 10 9 L Q Mpc -1 , about 
a factor of 2 larger th an the previous determinatio n at 
z ~ 2. Note that the IPerez-Gonzalez et al.l (|2005j ) de- 
termination of the IR LD at z ~ 2.4 is consistent with 
our measurement, primarily for two reasons. First, they 
only excluded the most extreme AGN, which by num- 
ber made up 5% of their sample, whereas other sur- 
veys indicate a larger contamination rat e of 10 — 20% 
(see § 3.2 of Perez-Gonzalez ct al. 2005 and references 
therein). Second. [Perez-Gonzalez e~a l. (2005) use a con- 
version relation between rest-frame 12 /im and IR lu- 
minosity that has been shown to overproduce by fac- 
tors of > 2 the IR luminosities of the most IR lumi- 
nous galaxies based o n 70 and 160 /im stacking analysis 
dPapovich et al.ll2007f) . This result is also supported by 
I Daddi et al.l (|2007ph . who also find a systematic excess 
of rest-frame 8 /jm emission relative to UV and radio 
emission for a sampl e of z ~ 2 star-forming ULIRGs. 
Further, I Daddi et al.l (|2007al ) suggest that the obscured 
AGN fraction among ULIRGs at z ~ 2 is significantly 
larger than previously inferred. A more conservative ap- 
proach of e xcluding IR-luminou s AGN, e.g., such as that 
adopted bv lCaputi et al.1 (|2007h . and assuming their cal- 
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Fig. 27. — IR luminosity den sity as a funct ion of redshift, including data from Yun ct al. (2001) at z = 0, Flor es et al.l 11199 9) at 
0.2 < z < l .O.IBareer et al.l ||200(1) at z ~ 2, and lChapman et al.l (I2005T ) at z > 2; all shown with open squares. Results from lCaputi et al.l 
(2007) and P crcz-Gonzalcz ct al. (2005) are shown by the open circles and gr een lines, respectively. Our points at 1.9 < z < 2.7 and 
2.7 < 2 < 3.4 are shown by the solid pentagons, assuming the Caouti ct al. ( 2007) calibration between i/Li/(8/im) and Lib,- ThelReddv et ah 
(2006b) calibration would raise the points by « 30%. The ULIRG, LIRG, and total contributions to the IR LD inferred bv lLe Floc'h et al 
(2005) are indicated by the light, medium, and dark gray shaded regions, respectively, up to z = 1. The inferred contribution of LIRGs 
from our analysis at z ~ 2 — 3 is denoted by the shaded rectangles. 



ibration between 8 fim and IR luminosity, results in an 
IR LD from U LIRGs that is roughly a fac tor of 2 — 3 
lower than the IPerez-Gonzalez et al.1 (|2005| ) determina- 
tion. Th e critical point of this discu ssion is that while 
both t he IPerez-Gonzalez et al.l (|2005l ) and ICaputi et al.1 
(|2007f l studies were able to probe the ULIRG regime 
of the IR LF with a high degree of completeness, the 
shallowness of their data precluded strong constraints on 
galaxies significantly fainter than 10 12 L Q . Our sample 
and analysis yield the first constraints on the IR and 
bolometric LFs for moderately luminous galaxies at red- 
shifts 1.9 < z < 3.4, thus allowing us to evaluate the 
unobscured SFRD over a larger range of intrinsic lumi- 
nosity and redshift than previously possible. Our results 
suggest that the luminosity density contributed by sub- 
ULIRGs with Lm < 10 12 L Q at z ~ 2 is larger than 
previously inferred. Table [10] indicates that such galax- 
ies would comprise roughly 70% of the total IR LD at 
z~2. 

It is important to take a step back at this point and re- 
examine the systematic effects induced by the unknown 
attenuation distribution of UV faint galaxies, since this 
is the dominant uncertainty in our determination of the 
faint end of the IR LF. The only way to reconcile the 
relatively steep faint-end slope of the UV LF (a = —1.6) 
with the shallow faint-end sl ope of the IR LF as sug- 
gested bv lCaputi et al.l (|2007l ) and lPerez-Gonzalez et al.1 
(2005) is if there is a sudden change in the attenuation 



properties of UV galaxies such that those with 1Z > 25.5 
have negligible dust attenuation. The implication is that 
such UV-faint galaxies would be forming stars in pristine 
metal-poor ISM with very low dust-to-gas ratios, akin 
to "dwarf" galaxies in the local universe. While there 
are certainly likely to be such dwarfs at high redshift, 
they would have to dominate the number counts on the 
faint-end of the UV LF. Future deep spectroscopic obser- 
vations should place both the determination of the UV 
faint-end slope and the extinction properties of sub-L* 
galaxies on a more secure footing. However, as we dis- 
cuss in the next section, comparison of our total IR LDs 
with those correspond ing to SFRD valu es estimated in 
previous studies (e.g. JReddv et afll2005h suggest the IR 
LD cannot be much lower than the value derived here, 
effectively placing a constraint on the attenuation of sub- 
ULIRGs. 

Lastly, we note that UV emission comprises a non- 
negligible fraction of the bolometric luminosity of LIRGs. 
This is an effect that becomes more pronounced at higher 
redshift as the average dust attenuation of galaxies of 
a given bolometric luminosity decreases, as already dis- 
cussed. The bolometric LFs derived in § [6] allow us to 
assess the contribution of galaxies to the total LD, tak- 
ing into account both the obscured (IR) and unobscured 
(UV) luminosity densities. The total bolometric lumi- 
nosity density can be calculated by integrating our bolo- 
metric LF from L ~ to L = 10 12 L Q , and adding the 
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contribution of ULIRGs from ICaputi et al.l (|2007t) . As- 
suming that the fraction of bolometric luminosity emer- 
gent in the UV in ULIRGs is negligible, then the to- 
tal bolometric LD at z ~ 2 is « 1.5 x 10 9 L Q Mpc -3 . 
Roughly 80% of this bolometric LD arises from galaxies 
with Lboi < 10 12 L Q . This bolometric LD is more than 
25% larger than what we would have inferred from the 
IR LF alone because the former includes the contribution 
of the LD emergent at UV wavelengths. 

To summarize, there are essentially four points worth 
keeping in mind from our analysis of the IR and bolomet- 
ric LD. First, the evolution of the UV LD shows a marked 
difference from the evolution of the dust-corrected Ha 
and IR LDs between z ~ 2 and z ~ 0. The difference 
can be explained by an evolution in the average extinc- 
tion correction between z ~ and z ~ 2. Second, we 
find that the IR LD at z ~ 2 is at least as large as the 
value at z ~ 3, implying that the decline in SFRD to 
the local value must have occurred after z ~ 2. Third, 
while our analysis becomes increasingly incomplete for 
the most luminous galaxies at z ~ 2 — 3, it does provide 
the first spectroscopic constraints on the moderate lu- 
minosity (e.g., LIRG) regime of the IR LF. Even taking 
into account the significant uncertainties associated with 
the dust obscuration of UV-faint galaxies, these results 
suggest that previous studies have significantly underes- 
timated the contribution of galaxies with Lir < 10 12 L Q 
to the IR luminosity density. Finally, taking into ac- 
count the emergent UV luminosity density of galaxies, 
we find that sub-ULIRG galaxies comprise roughly 80% 
of the total bolometric LD at z ~ 2. In the next section, 
we will discuss these results in the context of the global 
SFRD. 

8.5. Constraints on the Global Star Formation Rate 
Density 

We have converted the results of Figure [231 and [27] to 
star formation rates (SF Rs) using t he | Kennicuttl (1998) 
relations and assuming a ISalpeter] (|1955f ) IMF from 0.1 
to 100 Mq, with the results summarized in Figure l28l 25 
The UV points at low redshift (z < 1) are corrected 
assum ing a factor of 7 attenuation ( Schimin ovich et al.l 
|2005[) . Our UV points are correc t ed by an average factor 
of 4.5 (e.g., IReddv et al.ll2006bl: iReddy fc Steide1 [200l 
iNandra et all l2002t ISteidel et alj I1999D . Note the dif- 
ference in extinction correction between the low redshift 
points and our z ~ 2 — 3 determinations, again reflect- 
ing the dependence of extinction on redshift. The higher 
redshift points of Bouwens et al. (2007, submitted) are 
corrected according to the extinction factors published 
therein. While there are a number of systematics that 
can affect our interpretation of SFRD plots such as Fig- 
ure[28] (e.g., changes in the IMF with redshift), our anal- 
ysis has allowed us to carefully and quantitatively assess 
two of these systematics: accounting for sample incom- 
pleteness and extinction. Our measurements indicate an 
SFRD at z ~ 2 that is at least as large as the value at 
z ~ 3. Applyi ng a factor of 4.5, as suggested by stacking 
analyses (e.g. IReddv fc Steide1l2004 IReddv et aLll2006bt 
INandra et alj|2002h to correct our UV estimates for ex- 
tinction yields values that are in general accord with the 

5 Assuming the more realistic Chabrier ( 2003) IMF will reduce 
the SFRD by a factor of ~ 1.8. 



IR estimates. On the one hand, we might not have ex- 
pected such good agreement given that the factor of 4.5 
has only been measured for LIRGs, and the average fac- 
tor applied to the total UV LD may be lower depending 
on the extinction of UV-faint galaxies. However, we note 
that the IR-estimated SFRD includes the contribution 
from ULIRGs (whereas the UV-estimated SFRD does 
not explicitly take them into account), so the difference 
in the two estimates is less than we might have expected. 
The important point is that despite the significant uncer- 
tainties regarding the attenuation of UV-faint galaxies, 
applying a factor of 4.5 to total UV LDs as advocated in 
many studies is not too far off from the value obtained 
from IR estimates. 

As alluded to previously, independent estimates of the 
SFRD based on "ce nsus" studies of high redshift galaxies 
(|Reddv et al.| [2005) can be compared with our complete- 
ness corrected estimate. Such census studies estimate the 
global SFRD by adding up the contribution from galaxies 
dire ctly targeted b y var ious selection criteria. For exam- 
ple, IReddv et all (j2005[ ) estimated the SFRD at z ~ 2 
by taking into account the overlap between ga laxies se- 
lected using the BX BzK (iDaddi et al.ll2004bf l. Distant 
Red Galaxy (DRG; iFranx et al.ll20o¥ criteria, and sub- 
millimeter selection (e.g.. Ichapman et al.ll2005l ). Taking 
into account th e overlap between ga laxies selected using 
these methods, IReddv et al.l (|2005h compute an SFRD 
at z ~ 2 of 0.15 ± 0.03 M Q yr" 1 Mpc -3 for galaxies 
with K < 25.5, K S {AB) < 23.8, including the contribu- 
tion from submillimeter-bright (Ssso/x ~ ^ m Jy) galax- 
ies. Since the total SFRD cannot be smaller than that 
derived from census studies of directly detected (optical, 
near-IR, and submm bright) galaxies, it suggests that 
our SFRD measurement cannot be significantly lower 
than the value derived here, ~ 0.2 M Q yr _1 Mpc -3 . 
This in turn implies that the attenuation of sub- ULIRGs 
cannot be so low as to bring down our total SFRD 
estimate to the point where it is in violation of the 
census-computed SF RD. As an example, integrating the 
ICaputi et all (|2007f ) IR LF at z ~ 2 yields an IR LD of 
~ 6.6l^ x 10 8 L Mpc" 3 , corresponding to an SFRD of 

0.11±°; 2 M yr" 1 Mpc" 3 . This "total" value is already 
comparable to, if not smalle r, than the census value of 
0.15 ± 0.03 M Q yr" 1 Mpc" 3 (|Reddv et al.ll2005h : the lat- 
ter can be treated as a lower limit to the SFRD. The 
implications are that the faint-end of the IR LF is un- 
likely to be as shallow as that predicted from previous 
IR surveys, that such UV-faint galaxies are likely to have 
non-negligible dust attenuation, and that the contribu- 
tion of sub-ULIRGs to the total LD and SFRD must be 
larger than previously inferred. Finally, we note that 
there are no observational constraints on the presence 
of dusty LIRGs at z ~ 2 — 3 that are both faint in the 
UV and fall below the detection threshold of the Spitzer 
surveys. However, large numbers of such UV-faint dusty 
LIRGs would serve to increase the average attenuation of 
UV-faint galaxies and would strengthen our conclusions 
regarding the increased contribution of sub-ULIRGs to 
the total LD and SFRD at z ~ 2 - 3. 

One of the principle conclusions from this work is that 
galaxies with Lin ~ ibol ^ 10 12 L Q account for w 70% 
of the SFRD at 1.9 < z < 2.7. The fraction rises to 
~ 80% if we take into account the emergent UV luminos- 
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Fig. 28. — Summary of UV (blue open circles) and IR (red open squares) estimates of the global star formation rate density as a function 
of redshift. The UV points have been corrected for extinction in a differential manner (i.e., as a function of redshift): determinations at 
z < 1 are corrected by a factor of 7; our values at z ~ 2 and z ~ 3 are corrected by an average factor of 4.5; the high redshift (z > 4) 
determinations are corrected by the factors given in Bouwens et al. (2007, submitted): the dust corrections at z ~ 4 — 5 are determined 
from a linear interpolation of the z ~ 3 and z ~ 6 values. Our dust-corrected UV determinatio ns a re indicated by the large blue pentagons. 
Our IR estim ates, computed using a differential extinction recipe with UV magnitude (see § I6.3|l and which includes the contribution of 
ULIRGs from [Caputi et al. (2007) at z ~ 2, are indicated by the red open pentagons. The errors on our UV and IR determinations at 
z ~ 2 and z ~ 3 are ~ 20% (derived from the errors on the luminosity density) and are smaller than the size of the large pentagons. 



ity density of sub-ULIRGs. Our analysis suggests that 
much of the star formation activity at z ~ 2 may take 
place in faint and moderately- luminous galaxies. As a 
consequence, such faint to moderate luminosity galax- 
ies must have accounted for a significant fraction of the 
stellar mass density formed between z ~ 3 and z ~ 2. 
Assuming a constant SFRD of ~ 0.16 M© yr _1 Mpc -3 
(i.e., the average of our IR estimates of the SFRD at 
z ~ 2 and z ~ 3) between z = 3.4 and z — 1.9 (the 
limiting redshifts of our analysis) implies that ~ 45% 
of the total ste llar mass density of the Universe (e.g., 
ICole et aT1l200lD formed b etween these redshifts (e.g., see 
also lDickinson et al.ll2003l ) . If the fraction of the luminos- 
ity density contributed by galaxies with Lboi ^ 10 12 L Q is 
roughly constant (w 80%) between z = 3.05 and z = 2.3, 
as suggested by our analysis, then it implies that one- 
third of the present-day stellar mass density was formed 
in such galaxies. 26 

Finally, including the z > 4 rest-frame UV deter- 
minations (no such corresponding Ha or IR data ex- 
ist at corresponding redshifts), sug gests that the SFRD 
falls off at these early epochs (e.g.. iBouwens et aDl2006l 

26 We assume a Salpeter IMF in these calculations. A Chabricr 
IMF will reduce the stellar mass density estimates by a factor of 
~ 1.8, but the relative contribution of sub-ULIRGs between z = 
3.05 and z = 2.3 to the present-day stellar mass density will be the 
same, assuming the IMF does not evolve between z ~ 2 — 3 and 
the present-day. 



and references therein). Assuming a constant IMF im- 
plies that the SFRD at z ~ 2 is a factor of ~ 9 
larger than that observed at z ~ 6 (corrected for ex- 
tinction). This trend may be suggesting a hierarchical 
buildup of galaxies at early times. The important re- 
sult of our analysis is that this "growth" appears to halt 
and stay constant for roughly 1.2 Gyr, between z ~ 4 
and z ~ 2 (Figure [23]) . This result is not far from 
expectations if the rate of galaxies that are commenc- 
ing rapid star formation and populating the bright-end 
of the UV LF is offset by the rate of galaxies that are 
fading because of gas exhaustion or some other trunca- 
tion of the star formatio n such as through AGN or su- 
pernovac fee dback (e.g . . iKriek et alj 120061: iReddv et all 
I2006U 120051 : lErb et al.ll2006cf l . In fact, this very epoch 
hosts the emergence of a significant population of qui- 
escent galaxies (e.g.. iFranx et alj 120031: iRudnick et al.l 



20031: Ivan Dokkum et alj|2004 IReddv" et al. 2005, 2006b; 



Krick etal .1 120061 ). The balance between rapidly star 



forming galaxies and those that are fading appears to 
saturate around z ~ 1 — 2, at which point the latter be- 
comes the dominant effect, leading to a decrease in the 
number density of bright galaxies between z ~ 2 and 
z = (Figure [53]). As discussed above, this reversal 
in the evolution of the UV LF is likely due to gas ex- 
haustion resulting from any number of processe s (e.g., 
SN-driven outflows, AGN feedback, see also iBell et all 
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2005) . Once gas exhaustion has occurred and star forma- 
tion proceeds quiescently, at least a couple of mechanisms 
have been proposed to explain why gas is unable to cool 
onto galaxies wi th the largest stella r masses, including 
AGN feedback (|Croton et al.l 120061 : iScannapieco et all 
l2005t iGranato et al.ll2004h and dilution of infalling gas 
due to virial shock ( Dekel fc Bir nboim 2006h. It is inter- 
esting to note that it is around thi s epoch, z ~ 2. that 
AGN activity appears to peak (e.g . , iHopkins et al.ll2007t 
iFan et al.l[200lt iShaver et al.lfl996fl . 

9. CONCLUSIONS 

We have used the largest existing sample of spectro- 
scopic redshifts in the range 1.9 < z < 3.4 to evaluate 
the luminosity functions (LFs) of star-forming galaxies 
at rest-frame UV, Ha, and infrared (IR) wavelengths. 
The sample of rest-frame UV selected galaxies includes 
~ 15000 photometric candidates in 29 independent fields 
(for a total area of ~ 0.9 square degrees) of which ~ 2000 
have spectroscopic redshifts 1.9 < z < 3.4. The large 
spectroscopic database yields critical constraints on the 
contamination fraction of our sample from objects at 
lower redshifts (z < 1.4) and AGN/QSOs; statistics that 
are vital to accurately estimate the bright-end of the LFs. 
We use our extensive sample to correct for incomplete- 
ness and recover the intrinsic rest-frame UV LF at z ~ 2 
and z ~ 3. Combining this result with Ha and Spitzer 
MIPS data in several of our fields enables us to infer the 
Ha and IR LFs of star-forming galaxies. The principle 
conclusions of this work are as follows: 

1. The fraction of star- forming galaxies with rest- 
frame Lya equivalent widths Wiyc > 20 A in emission 
(/20) increases with redshift in the range 1.9 < z < 3.4, 
independent of selection bias, from a value of 8% at 
1.90 < z < 2.17 to 23% for Lyman-break galaxies at 
2.70 < z < 3.40. If the general expectation is that young 
and less dusty galaxies show Lya in emission, then the 
trend of increasing /20 with redshift reflects the decrease 
in average galaxy age and metallicity with increasing red- 
shift. 

2. Based on integrating our maximum- likelihood LFs, 
the fraction of star-forming galaxies with redshifts 1.9 < 
z < 2.7 and M A b(1700A) < -19.33 that have colors that 
satisfy BX criteria is » 58%. Similarly, the fraction of 
star-forming galaxies with redshifts 2.7 < z < 3.4 and 
M A b(1700A) < -20.02 that have colors that satisfy the 
LBG criteria is « 47%. The total fraction of 1.9 < z < 
3.4 galaxies with 1Z < 25.5 that satisfy either the BX 
or LBG criteria is 55%. We find little evolution in the 
rest- frame UV LF between z ~ 3 and z ~ 2. Correcting 
for extinction implies the dust-corrected UV luminosity 
density at z ~ 2 is at least as large as the value at z ~ 
3, and roughly 9 times larger than the value at z ~ 6. 
This evolution reverses at redshifts z < 2, where gas 
exhaustion is likely to dominate the evolution of UV- 
bright galaxies. 

3. The incompleteness-corrected estimates of the 
E(B — V) distribution indicate very little evolution in 
the average dust extinction of galaxies between z ~ 3 
and 2^2, and that such a distribution is approximately 
constant to our spectroscopic limit of 1Z = 25.5. These 
results are in agre ement with stacked X-ray analyses of 
z ~ 2 - 3 galaxies (|Reddv k. Steidell 12004 iNandra et all 
2002). However, examining the attenuation distribution 



of galaxies over a larger dynamic range in lookback time 
indicates an increasing extinction per unit star forma- 
tion rate with d ecreasing redshift (jReddv et alj l2006bl : 
iBuat et "aill2007h . The implication of this trend in aver- 
age extinction is that the evolution in the dust-corrected 
UV LD is less and more pronounced, respectively, at low 
(z < 2) and high [z > 2) redshift than what one would 
have predicted from the observed UV LD (e.g., see also 
iBouwens et~al]|2006l ). 

4. Factoring in the contamination rate of our sample 
from galaxies at lower redshifts and AGN/QSOs with 
redshifts 2.7 < z < 3.4, we find no evidence for an excess 
of UV-bright galaxi es over what was inf erred in the ini- 
tial L BG studies of ISteidel jtap (| 19991) and iDickinsonl 
(1998), as has been recently claimed. 

5. The incompleteness-corrected rest- frame UV se- 
lected sample and deep Spitzer MIPS data in multi- 
ple fields are combined to yield the first spectroscopic 
constraints on the faint and moderate luminosity sub- 
ULIRG (Ltr < 10 12 L Q ) regime of the total infrared 
luminosity functions at z ~ 2 and z ~ 3. We use 
this information to show that the number density of 
-^IR ^ 10 12 L© galaxies has been significantly under- 
estimated by previous studies that have relied on shal- 
lower IR data. After accounting for the emergent UV 
luminosity, and assuming a realistic range of attenua- 
tion for UV faint galaxies, we find that w 80% of the 
bolometric (IR+UV) luminosity density at z ~ 2 comes 
from galaxies with Lboi < 10 12 L Q . Assuming a con- 
stant SFRD of 0.16 M yr" 1 Mpc -3 between z = 3.4 
and z = 1.9 (the limiting redshifts of our z ~ 2 — 3 
analysis), suggests that Lboi < 10 12 L Q galaxies at these 
redshifts were responsible for approximately one-third of 
the buildup of the present-day stellar mass density. 

6. Our estimate of the total SFRD at z ~ 2 is consis- 
tent with the lower limit on the SFRD provided by cen- 
sus studies of optical, near-IR , and submill i meter bright 
galaxies at similar redshifts (|Reddv et al.l [2005) . Our 
analysis takes into account the systematics associated 
with the attenuation of UV-faint (1Z > 25.5) galaxies. 
Assuming more extreme changes in the dust attenuation 
of UV-faint galaxies than considered here would be re- 
quired to reconcile the steep and shallow faint-end slopes 
of the UV and IR LFs, respectively, implying that the 
vast majority of UV-faint galaxies would be forming stars 
from chemically pristine gas. However, such an extreme 
scenario would result in a total SFRD that is compara- 
ble to, i f not smaller than , the lower limit from census 
studies (jReddv et al.ll2005T ). The implications are that 
some significant fraction of sub-ULIRGs must have non- 
negligible dust extinction, the faint-end of the IR LF 
must be steeper than what previous studies have sug- 
gested, and the contribution of sub-ULIRGs to the total 
SFRD must be significantly larger than previously in- 
ferred. 
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